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ABSTRACT 

A  five  degree-of-freedom  linear  model  for  the  XR-3 
surface  effect  ship  is  developed  for  constant  speed  operation, 
The  weight  removal  transient  response  and  the  roll  behavior 
of  the  craft  in  turn  maneuvers  are  investigated  by  both 
linear  and  six  degree-of-freedom  nonlinear  models,  and  the 
linear  model  simulation  results  are  compared  with  the  non- 
linear model  simulations.   Some  of  the  parameters  of  the 
craft  are  investigated  for  roll  motion. 
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I.   INTRODUCTION 

A.  BACKGROUND 

The  Surface  Effect  hybrid  type  high  performance  ship  has 
been  of  great  interest  recently,  and  in  1970  the  XR-3  test 
craft  was  delivered  to  the  Naval  Postgraduate  School.   As 
seen  in  [7]  Leo  and  Boncal  converted  the  digital  simulation 
program  for  the  100-B  test  craft  to  represent  the  XR-3,  thus 
providing  a  nonlinear  model.   In  [1],  Gerba  and  Thaler  have 
developed  a  heave-only  linear  model  and  further  studies  were 
made  by  Barnes  in  [2]  by  developing  a  two  degree-of -freedom 
linear  pitch  and  heave  model  representing  the  XR-3  test  craft. 

B.  OBJECTIVE 

In  order  to  investigate  the  roll  behavior  of  the  XR-3 
test  craft  during  turn  maneuvers,  a  five  degree-of -freedom 
linear  model  is  developed  in  this  thesis.   Some  modifications 
are  made  for  pitch  and  heave  modes,  therefore  test  of  the 
vertical  plane  by  means  of  weight  removal  is  reinvestigated 
and  compared  with  the  nonlinear  simulation  results. 

In  [3] ,  the  roll  behavior  of  the  XR-3  and  the  sidewall 
sensitivity  are  studied  using  a  six  degree-of -freedom  non- 
linear model.   Similar  studies  are  made  in  this  thesis  after 
developing  a  five  degree-of-freedom  linear  model  and  comparing 
with  the  nonlinear  model  results.   Finally,  some  modifications 
are  suggested  for  nonlinear  modeling. 
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II.   NONLINEAR  XR-3  MODELING  FOR  FIVE  DEGREE-OF-FREEDOM 

A.   ASSUMPTIONS 

The  following  assumptions  and  simplifications  are  used 
during  development  of  a  constant  speed  XR-3  test  craft  model: 

1.  Aerodynamic  drag  on  the  forward  seal,  and  superstruc- 
ture skin  friction  forces  are  neglected. 

2.  Momentum  (ram)  drag  is  disregarded  for  all  modes. 

In  other  words,  it  is  accepted  that  the  location  of  the  lift 
fans  with  respect  to  the  craft  C.G.  could  not  affect  signi- 
ficantly the  roll  and  pitch  stability.   Also,  air  stagnation 
force  in  the  z-direction  is  assumed  to  be  negligible. 

3.  During  dynamic  response  of  the  craft,  only  air 
leakage  with  constant  leakage  area  at  the  stern  seal  is 
assumed.   No  other  air  escape  is  considered  such  as  possible 
leakage  from  sidewalls  or  bow  seal. 

4.  Plenum  pressure  and  air  mass  thermodynamics  are 
represented  by  ideal  diabatic  processes  with  the  air  constant 
specific  heat  ratio  y  =  1.4. 

5.  Effective  plenum  roof  area  is  developed  as  a  function 
of  craft  weight. 

6.  Planing  forces  are  concentrated  at  the  centroid  of 
the  keel  longitudinal  cross  section  through  the  fluid  medium. 

7.  Pitch  and  roll  damping  moments  due  to  added  mass  are 
included  in  the  craft  dynamics,  but  their  contribution  to 
heave  motion  is  neglected. 
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8.  The  thrust  forces  and  all  hydrodynamic  resistive 
forces  (drag)  in  the  z-direction  cancel  each  other,  since 

at  constant  speed  operation,  these  forces  are  in  equilibrium. 
However  they  might  contribute  to  the  moment  equations  such 
as  pitch  moment,  but  considering  their  effective  moment  arms 
with  respect  to  the  C.G.  are  almost  equal  to  the  thrust  moment 
arm,  so  that  their  contribution  in  the  moment  equations  is 
also  neglected. 

9.  For  seal  modeling,  vertical  resistive  forces  due  to 
hydrodynamic  response  are  taken  into  account  including 
corresponding  effective  wetted  seal  area. 

10.   Simplified  craft  geometry  is  used  for  plenum, 
buoyancy  and  seal  modeling. 

B.   COORDINATE  SYSTEM  AND  SIGN  CONVENTIONS 

1.  The  cartesian  coordinate  system  is  used  with  its 
origin  located  on  the  water  line  corresponding  to  the  craft 
center  of  gravity. 

2.  The  x-axis  lies  on  the  water  line  with  the  positive 
longitudinal  displacement  being  measured  forward. 

3.  The  y-axis  lies  in  the  vertical  cross  section  plane 
with  the  lateral  dislacement  being  measured  positive  to 
starboard. 

4.  The  z-axis  lies  in  the  vertical  cross  section  plane 
with  the  positive  vertical  displacement  being  measured 
downward. 


13 


5.  z  represents  heave  C.G.  displacement  measured  from 
the  origin   (at  water  line)  to  the  center  of  gravity. 

6.  z   is  one  of  the  craft  dimensional  values  measured 

s 

from  the  keel  to  the  C.G. 

7.  Pitch  is  the  rotational  motion  about  the  y-axis, 
angular  pitch  displacement  is  represented  by  9 ,  positive 
pitch  angle  is  measured  upward  (boat  noses  up) . 

8.  Roll  is  the  rotational  motion  about  the  x-axis, 

cj)  angle  represents  angular  roll  displacement,  positive  roll 
angle  is  measured  to  starboard  (boat  heels  to  starboard) . 

9.  Similarly  the  yaw  is  the  rotational  motion  about  the 
z-axis,  positive  yaw  is  measured  in  the  clockwise  direction. 

10.  Positive  roll,  pitch  and  yaw  moments  are  to  be 
defined  as  the  moments  in  the  positive  roll,  pitch  and 
yaw  directions  correspondingly. 

11.  As  a  reference  zero  pitch  and  roll  angle  provide  the 
x-y  plane  parallel  to  the  water  surface. 

The  craft  coordinate  system  is  shown  in  figure  1. 

C.   SIMPLIFIED  CRAFT  MODELING  AND  GEOMETRY 

1.   General  Simplified  Geometry  and  Dimensions  of  the 
XR-3  Test  Craft 

In  figure  1  the  sideview  simplified  geometry  of  the 

XR-3  test  craft  is  shown.   Also  in  figure  3  the  top  view  is 

illustrated.   The  variables  and  dimensions  in  these  figures 

are  defined  as  follows; 
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Bow  Seal 
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Figure  2.   Bow  seal  geometry 
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Figure  3.   Simplified  geometry  (top  view) 
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I ,   =  Draft  of  the  craft  at  center  of  gravity 

C   =  Plenum  pressure  center  (the  centroid  of  the 
"     effective  plenum  area) 

£_,£,  =  The  longitudinal  distance  of  the  sections' 
centroids  from  the  C.G.  (forward  and  aft 
section  centroids  respectively) 

£_,£.   =  Bow  and  stern  seals  front  face  distances 
at  water  level 

I-      =   Centroid  distance  of  the  planing  force 
effective  area 

Z^,i2     =  Forward  and  aft  sections  sidewall  submerged 
region  length. 

All  craft  dimensional  values  are  taken  from  the  XR-3 
six  degree-of-freedom  model  scale  drawings  as  indicated  in 
[2]  and  [3]. 

2.   Sidewalls 

As  reported  in  [2] ,  the  sidewall  modeling  is  developed 
from  the  SIDEWALL  subroutine  of  the  6  DOF  computer  model. 
Like  craft  dimensions,  all  sidewall  dimensions  and  their 
related  geometric  configurations  are  based  on  scale  model 
drawings.   For  simplicity  of  calculation,  the  sidewalls  are 
divided  into  two  sections  called  forward  and  aft.   During 
dynamic  response  of  the  sidewalls,  these  sections  are  con- 
sidered individually  with  their  own  centroids,  so  related 
calculations  can  be  done  easily.   Especially,  the  sidewall 
modeling  is  developed  in  order  to  determine  buoyancy  forces 
and  moments.   By  means  of  this  modeling,  the  average  sidewall 
width,  average  draft  at  section  centroid  and  effective 


17 


plenum  area  are  also  determined  and  used  extensively  in  the 
rest  of  the  calculations. 

The  sidewall  vertical  cross  sections  at  the  section 
centroid  forward  and  aft  are  sketched  in  figures  4(a)  and  4(b) 
respectively.   From  these  figures,  the  average  forward  side- 
wall  width  at  the  buoyancy  center  is  found  to  be 


Ff  =   2  tan  B,  +  bf  (II-1) 


where 


I-     =  Draft  at  forward  buoyancy  center 

tan  Bj   =  Tangent  of  the  sidewall  deadrise  angle 
at  the  forward  buoyancy  center 

bf  =   Sidewall  flat  surface  width  at  the  kell 
(forward  section) 


Therefore  submerged  sidewall  volume  of  the  forward 
section  is  given  as  follows 


Vl  =  ^f  lf   ll  (H-2) 


£f  can  be  formulated  as 


lc     =  I-   -  l-   tan  8 
r       a     5 


then  the  small  angle  approximation  tan  6=0  gives 
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Figure  5.   Plenum  modeling 
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li  ld   ~  £5 


(II-3) 


As  defined  before,  £,  draft  at  C.G.  is 


%,      =   z   +  z 
a      s 


(H-4) 


Similar  expressions  can  be  derived  for  the  aft 


section, 


2  tan 


+  b 


(II-5) 


where 


V, 


sidewall  flat  surface  width  at  the  keel 
(aft  section) 


F   A.  I 
a  2   a 


£.  +  lc    0 
a     6 


(II-6) 


Note  that  &-,    i_   are  functions  of  draft  and  pitch 
angle. 

3.   Plenum  Modeling 

In  [2]  plenum  modeling  is  presented  including  plenum 
pressure  variation  due  to  changing  of  draft.   In  the  develop- 
ment of  a  pitch,  heave  and  roll  linear  model,  the  simplified 
plenum  geometry  and  dynamics  which  were  derived  by  Barnes 
are  considered  and  used  in  our  plenum  dynamics  representation. 

The  variables  and  dimensions  for  plenum  modeling, 
as  seen  in  figure  5,  are  found  as 
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£P   ■   £s  "  S  (Tr}  (II"8) 


*b    ■     bw  lP  (II"9> 


*9      £8  "  Zo 


where 


I        =  Length  of  plenum  at  the  water  line 

A,   =   Effective  plenum  area 

b   =  Lateral  distance  from  one  sidewall  to 
another  in  the  plenum  region. 


£  ,  £„,  £,,  £   and  h  are  all  constants  related  to 

craft  dimensions.   Only  the  effective  plenum  area  A,  varies 

when  the  craft  draft  varies,  that  means  variation  of  the  plenum 

pressure  centroid  C   is  forward  or  backward. 

P 

4.   Bow  and  Aft  Seals 

The  seals  are  an  involved  additional  dynamic  system 
with  their  own  degree  of  freedom  relative  to  the  craft,  and 
they  carry  little  load  compared  with  the  other  components 
of  the  craft.   However  in  the  simplified  model,  the  seal 
motions  relative  to  the  craft  are  neglected. 

The  assumed  seal  geometry  is  shown  in  figure  2 , 
developed  in  [2] .   Using  an  appropriate  geometric  configura- 
tion, the  seal  wetted  area  is  derived  as  a  function  of  craft 
draft  and  pitch  angle. 
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From  [2] ,    the  width  of  the  wetted  surface  for  the 
bow  seal  is; 


1,-1-    tan  9 
xf  =     sin  31° (II-1°) 


where 


sin  31°   =   sin  of  the  angle  between  water  surface 
and  the  bow  seal  front  face. 


The  numerator  of  this  ratio  represents  draft  at  the 
bow  seal  front  face. 

A  similar  expression  for  the  stern  seal  is; 


%,   +1.    tan 
d    4 
x 


a        sin  32° 

Multiplying  x^  and  x  by  fa  ,  bow  and  stern  seals 
wetted  areas  are  obtained  as  follows; 


A,-   =   fa   x. 
f      w  f 


A   =   fa   x 
a      w  a 


5.   Rudder  and  Propulsion  Hydrodynamics 

The  rudder  hydrodynamics  is  calculated  from  [4] , 
inserting  rudder  geometry  dimensional  values  into  the  fin 
equations.   Actually  we  are  interested  in  the  lateral  force 
(sway)  on  the  rudders  for  a  given  rudder  order  deflection. 
This  dynamic  equation  is  given  in  [4]  as  follows; 
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~  a  d      v 

Y        =      p    u      A    (l+v~)    C    [6     -(l+-£)    -"-]  (Il-ll(a)) 

r  rhrrhu 

r  r 


where 


A   =   Rudder  platform  area  (for  one  rudder  only) 

d   =  draft  at  the  stern  =   z  +  z   -  x   8 
s  s    r 

h   =  Vertical  depth  to  bottom  tip  of  the  rudder 

C   =   Rudder  lift  coefficient 

2tt  R 

=  = — t— -o-   (for  fully  wetted  flow  regime) 
Ra  +  3 

R   =  The  rudder  aspect  ratio 

cl 

5    =   Rudder  deflection  (the  sign  convention  of  the 
rudder  is  to  be  positive  for  port  turn) 

V,   =  V  +  x  r-z   p 

h  r      s  c 

x   =   longitudinal  coordinate  of  rudder  geometric 
center  (always  negative  value) 

V   =   lateral  speed  (ft/sec) 

r   =  yaw  rate  (rad/sec) 

In  order  to  simplify  the  rudder  dynamic  equation, 
the  following  coefficients  are  defined  as 


d 
e         hr 


Ct  =   p  u2  Ar(l  +  fZ)    cr 
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Then  equation  (11-11 (a))  becomes 


Yr  =   Ct(6r  "  Ce  "£>  (U-ll(b)) 


Let  us  define  z   to  be  the  vertical  coordinate  of 

r 

the  rudder  from  the  C.G.  to  its  geometric  center,  the  rudder 
roll  and  yaw  moments  can  be  calculated  from 


Kr  =   -  Yr  zr  (.11-11  (O) 


Nr  =  Yr  Xr  (H-ll(d)) 


During  turn  maneuvers,  the  propulsion  thrust  vector 
is  also  deflected  according  to  the  given  rudder  order.   There- 
fore the  lateral  component  of  the  thrust  contributes  yaw  and 
roll  moments  as  well  as  a  sway  mode,  but  this  diturbance  must 
be  considered  as  a  step  input  to  the  linear  modes  with  constant 
magnitude,  therefore  the  total  propulsion  sway  force  is 


Y     =   2  X.  sin  6  (11-12) 

pr        t      r 


where  X  is  the  propulsion  thrust  corresponding  to  craft  speed. 

Similarly  the  roll  and  yaw  moments  generated  by  the 
propulsion  system  are 


Kpr  -  Ypr  Zp  CH-13(a)) 
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N     =   Y    x  (11-13  (b)) 

pr      pr  p  v  '  ' 


where  x  and  z   are  the  longitudinal  and  vertical  locations 
P      P 

of  the  propeller  respectively. 

D.   DYNAMIC  FORCES  FOR  FIVE  DEGREE-OF-FREEDOM  MODES 
1.   Seal  Forces 

Considering  the  seal  configuration  and  dynamics ,  the 
pressure  difference  across  the  wetted  portion  of  the  seal  is 


P,   =  P.  -  P 
b      b    a 


where 


—  2 

P.   =  Plenum  gauge  pressure  (lb/ft  ) 


P   =  Atmospheric  pressure 


Since  the  seal  is  in  equilibrium,  the  vertical  hydro- 
dynamic  force  on  the  seal  is 


Zsf  =  "  Pb  Af 


I.    -    £-.9 
Pb  *   V  sin  31° )  ^-^ 


For  the  stern  seal  hydrodynamics,  a  similar  expression 
can  be  found  as 


Zsa  =  "  (Pb  +  V  VlEr^  (II"15) 
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where 


P.   =  Additional  pressure  supplied  into  stern  air 
bags  in  order  to  reduce  leakage. 

2.  Buoyant  Forces 

The  forward  section  buoyant  force  is  determined  from 
equation  (II-2) ,  it  is 

zbf  =  "  2  p  g  vx 

here  p  is  the  density  of  the  seawater  and  g  is  the  gravitational 
constant. 

Finally,  the  forward  and  aft  section  forces  are' 

Zbf  =  "  2  p  9(2  tan  B,  +  fafUf  *1         (II"16) 

Zba  *  "  2  p  ^2  tan  32  +  ba)£a  l2  (II"17) 

3.  Planing  Forces 

In  addition  to  the  pitch  angle  6 ,  the  sidewall  flat 
area  at  the  keel  determines  the  planing  effect  of  the  craft. 
When  the  craft  pitches  up,  the  frontal  area  (through  the  fluid 
medium)  causes  a  lift-drag  force  in  the  z-direction. 

The  planing  force,  from  Reidel's  work  13],  is; 


1      A   2   •   a 
Z  .      =   -yTrpAu   sin  9 
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where 

A  =   flat  surface  area  at  the  keel 
u  =   speed  of  the  craft 

Using  the  small  angle  approximation  and  including  the 
loss  coefficient  C  ■,  suggested  by  Barnes  in  [2]  ,  the  planing 
force  is  redetermined  as 


Zpl  =   -  Cpl  tt  p  A  u2  0  (H-18) 


A  similar  expression  is  derived  for  sidewall  dynamics 
in  [4] ,  but  in  our  modelling  equation  (11-18)  is  used. 
4 .   Pitch  Damping  Forces 

The  sidewalls  of  the  surface  effect  ship  influence 
moderately  heave,  pitch  and  roll  modes,  so  in  [4]  the  hydro- 
dynamic  and  hydrostatic  forces  and  moments  acting  on  the  side- 
walls  of  uniform  triangular  cross  section  are  derived  with 
particular  application  to  vertical  plane  motion  and  lateral 
motion. 

Considering  constant  speed  operation  of  the  craft, 
neglecting  some  terms  in  the  equations  provided  from  [4] , 
the  vertical  force  acting  on  only  one  sidewall  is  found  as; 


ZsWl  *  A33  Xs  U  <*  t11"19) 
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where 


A.,.*   =  vertical  added  mass  calculated  at  the  stern 

=   (tt/8)  p  Bs2 

B    =  value  of  sidewall  waterline  at  the  sidewall 
stern 

x    =   stern  distance  from  C.G.  (negative  value) 
s 

q    =  pitch  rate  (rad/sec) . 

For  the  pitch  damping  moment,  the  same  equation  is 
used  by  Barnes  in  [2] . 

5.   Roll  Damping  Forces 

The  XR-3  test  craft  is  dynamically  similar  to  100-B 
surface  effect  ship.   Therefore  we  can  use  dynamic  equations 
calculated  in  [4]  for  the  XR-3  test  craft.   The  heave-roll 
damping  force  for  one  sidewall  can  be  written  as; 


Zsw2   "•  "  A33  U  bo  P  CII"20) 


In  the  above  equation,  b  is  half  of  the  craft  width 
and  it  must  be  used  (+)  for  starboard  sidewall,  (-)  for  port 
sidewall. 

Similarly  the  lateral  damping  force  due  to  sidewall 
hydrodynamics  is 


YsW;L  =   -  A22  U(v  +  xsr  "  2sp)         CII-21) 
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here  A22  is  the  lateral  added  mass  calculated  at  the  stern, 
p  is  the  roll  rate.   A«2  is  calculated  from  [4] ,  that  is 


A22   =   tt/5  p  ds2 


6.   Plenum  Pressure  Lift  Force 


From  [2]  and  [4] ,  the  plenum  pressure  lift  force  is 
determined  assuming  uniform  pressure  distribution.  This  is 
formulated  as 

Zp  -  -  *b  ¥b 


from  equations  (II-8)  and  (II-9)  follow 


I      I. 

Z    =   -  fa   P.  U_ ^rr-^)  (11-22) 

p        w  b   8     h 


E.   DYNAMIC  PLENUM  AIR  MASS,  PRESSURE  EQUATIONS 

Assuming  constant  aft  seal  leakage  area,  the  leakage  flow 
rate  (Qout)  is  given  by  Gerba  and  Thaler  for  the  XR-3  heave 


only  model  in  [1] ,  then  we  have 

2  P 


Qout  -   CnA£(r-^)V2  (II"24) 

a 


where 


C   =   Nozzle  flow  coefficinet 
n 


A.   =   Leakage  area 
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p        =     Atmospheric  air  density. 

EL 

The  linearized  air  flow  (Q.  )  into  the  plenum  is  pre- 
sented in  [5]  according  to  the  fan  properties  and  calculated 
in  [1]  as  follows, 


Q.    =   n(Q   -  k   P.  )  (11-24) 

in        o    p  b 


where  n  is  the  number  of  fans  supplying  the  air  to  the  plenum, 
Q  is  the  steady  state  air  flow  rate  corresponding  to  steady 
gauge  pressure,  k  is  (1/slope)  of  the  P,  vs.  Q.   curve  at 
the  equilibrium  point. 

From  equations (11-23)  and  (11-24),  the  plenum  air  flow 
rate  can  be  written  as 


.      dm, 

M.   =  -=-£  =   P  (Q-   -  Q  .J  (11-25) 

d  dt      a  vin   wout 


The  equation  (11-25  J  reveals  thatif  the  fan  flow  rate 
(Q.  )  equals  to  the  leakage  flow  rate  (Q   .),  then  the  plenum 
air  mass  (M.)  is  constant. 

The  absolute  plenum  pressure  is  calculated  from 


pb  -  VvTlH  (II"26) 

b  Ka 


In  the  above  equation,  V,  is  defined  as  plenum  volume 
and  its  simplified  form  is  calculated  from  [1] , 
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V 


b  =   (Vn  -  *b  V  (II"27) 


where  V  is  to  be  the  empty  plenum  volume. 

F.   NET  HEAVE  FORCE 

The  resultant  heave  force  is  the  summation  of  the  all 
external  dynamic  forces,  it  is 


Z      ■      Z    c   +   Z         +    Z,  ,    +    Z,       +Z.+Z         +Z 
sf  sa  bf  ba  pi  sw  p 


where 


Z  =      Z  +   Z 

sw  sw.  sw_ 

Z    also  includes  a  pair  of  heave  buoyancy  forces  due  to 

bW 

roll,  but  they  are  in  opposite  directions.   Therefore  these 
stabilizing  forces  are  only  taken  into  account  in  the  roll 
motion. 

G.   PITCH  AND  ROLL  MOMENTS 
1.   Seal  Moments 
a.   Pitch  Mode 

The  seal  force  is  modeled  by  Barnes  in  [2]  acting 
on  the  center  of  the  seal  wetted  length.  The  effective  pitch 
moment  arms  are; 

Jl  -  =   £_  -  x ,  Z,    (for  forward  seal)  (11-28) 


'sf     £3    Xf/2 


£    =  %.    +  x  /2    (for  aft  seal)  (11-29) 

Sa         4       cL 
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From  this  follows  the  forward  and  aft  seal  moments 


Msf  =   "  Zsf  *sf  (II"30> 


Msa  '     "  Zsa  £sa  (II"31> 


where  Z  ~   and  Z   have  already  been  defined. 

Si  Scl 

b.   Roll  Mode 

The  flexible  seals  can  contribute  little  to  roll 
stability.   Actually  the  seal  configurations  and  hydrodynamics 
are  relatively  more  complex.   In  [4] ,  the  seal  roll  moments 
are  calculated  by  combining  the  roll  contribution  of  all  seal 
elements. 

Considering  the  seal  response,  the  submerged 
region  of  the  seal  can  be  calculated  and  its  roll  moment  con- 
tribution can  be  formulated  in  the  following  equations; 


K-  =   -  b  3/12  p  g  x-  tan  <j>  (11-30) 


Ksa  =   "  bw3/12  p  g  xa  tan  *  (H-31) 


For  this  modeling,  the  roll  moment  arm  is  considered 
as  (2/3  fa  ) ,  that  is  the  centroid  of  the  triangle  cross  section 
of  the  submerged  portion. 
2.   Buoyant  Moments 
a.   Pitch  Mode 

The  pitch  buoyancy  moments  are  calculated  from 
equation  (11-16)  and  (11-17)  assuming  constant  moment  arms 
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£_  and  &,.   The  pitch  buoyancy  moments  are 

"bf   ■   "  Zbf  *5  (II"32) 

"ba  =   Zba  l5  (II"33) 

b.   Roll  Mode 

As  discussed  in  the  heave  mode,  the  effect  of  the 
sidewall  buoyancy  due  to  the  roll  is  to  be  a  stabilizing  moment. 
From  figures  4 (a)  and  4 (b)  the  buoyancy  volume  change  can 
be  approximated  as 

fa  +  fa       I  I 

AVf   *   (JL2—+  2  tan  ^   <ti5^  *  bf }  ^       (II"34) 

Finally,  roll  moment  is 


Kbf   =   -  2  p  g  AVf  (H-35) 


Similarly,  the  aft  section  roll  moment  is 


Kba   =   -  2  p  g  AVa 


3.   Planing  Pitch  Moment 

The  planing  effect  due  to  the  sidewall  flat  area  at 
the  keel  contributes  only  pitch  motion  as  follows; 


Mi   =   Z  .  £_ 
pi      pi   7 
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4.   Pitch  Damping  Moment 

This  is  calculated  from  equation  (11-19)  as 


2 

M     =   -  2  Z     x    =-2A__ux 
sw  sw.   s  33     s 


In  the  above  equation,  the  factor  "2"  appears  because 
there  are  two  sidewall  sections. 
5.   Roll  Damping  Moment 

Combining  equation  (11-20)  and  (11-21) ,  the  roll 
damping  component  is 


K    =   2  Z    fa   -  Y    z+K        (11-37) 
sw        sw.   o    sw.   s    wave 


where  K     represents  the  additional  roll  damping  term 
wave 

because  of  vertical  wave  generation  in  the  roll. 
From  [3]  and  [4],  K     is 

Wave 


K      =   -  32/tt  b  2  B   P  (11-38) 

wave  o   c 


Here,  B  is  the  sidewall  vertical  added  mass  given  as 


2      2 

B   =   tt/8  p  I   I,     Cotan  B  . 
c  d 


6.   Plenum  Pitch  Moment 

From  equation  (11-22)  and  defining  X   to  be  the 

cp 

plenum  pressure  centroid,  the  plenum  pitch  moment  is  given 
as 
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M    =   -  Z   X     =   A,  (P,   -  P  )  X 
p        p   cp      b  b    a   cp 


where  X   is  calculated  as  a  function  of  craft  draft  in  [2] 
cp 

as  follows; 


X    =   X   +  (I      £,)/2h 
cp      c     g   d  ' 


X   =   £«/2  +  X 
c      8'      e.g. 


7.   Plenum  Roll  Moment 

The  destabilizing  static  plenum  roll  moment  is 
generated  when  the  craft  rolls  to  one  side.  Since  the  effec- 
tive plenum  pressure  centroid  shifts  laterally  to  the  other 
side  tending  to  increase  rolling  of  the  craft.   The  plenum 
roll  moment  is 


K   =  A.  (P.  -  P  )  (z   -  £.)* 
p      b  b    as    d  Y 


8.   Rudder  and  Propulsion  Roll  Moments 

These  roll  moments  have  already  been  calculated 
from  equation  (11-11  (c) )  and  (11-13  (a) ) . 

H.   SWAY  FORCES 

1.   Cross  Flow  Drag  (Sidewall) 

During  the  turn  maneuver,  the  cross  flow  drag  is 
generated  acting  on  the  outboard  sidewall.   The  direction 
of  this  force  is  negative  while  turning  to  the  left. 
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Neglecting  the  lateral  drag  component  inside  of  the  sidewall, 
the  net  cross  flow  drag  can  be  calculated  from 


Dsw  =   "  7   P  si^^r]    Cd  Vr2  i{id   ±    bo^  (II"39) 


where 

V   =  v  +  xe  r  -  (ze  -  A,/2)p 
r  s       s    a 

C,   =  cross  flow  drag  coefficient 

Considering  the  outboard  surface  of  the  submerged 
sidewall  sections  (dead  rise  angle  effect), the  lift  force 
in  the  z-direction  can  be  produced  due  to  the  lateral  drag 
deadrise  projection.   In  [4] ,  the  importance  of  the  sidewall 
deadrise  angle  is  mentioned,  however  at  the  hydrodynamic 
section  in  this  reference,  the  lift  component  of  the  drag 
force  has  not  been  derived. 

In  this  modeling,  the  lift  force  (Z,)  and  lateral 
drag  component  is  calculated  assuming  both  have  the  same  drag 
coefficients  (C,  =  1.28).   They  are 


Y,   =   DeT7  sin  6  (11-40) 

a      sw 


Zd  =   Dsw  COS  B  (H-41) 
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From  equation  (11-40)  and  (11-21) ,  the  total  hydro- 
dynamic  sway  force  of  the  sidewall  is 


Y     =   (D    +  Y    )  sinB+  Y  (11-42) 

sw       sw    sw  '     p    sw.  ' 


2.   Centrifugal  Force 

The  centrifugal  force  acting  on  the  C.G.  of  the  craft 
is  always  in  the  outboard  direction  during  the  turn  maneuver. 
Y  is  the  centrifugal  force,  given  as 

yc   -   "mir 

R,  the  turning  radius  is  equal  to  (u/r) ,  combining  this 
with  the  above  equation  gives 


Y   =   -  m  u  r  (11-43) 

c 


a  negative  sign  is  used  due  to  the  coordinate  system  sign 
convention.   Since  the  right  turn  produces  (+)  yaw  rate,  then 
Y  is  calculated  as  (-) ,  so  outboard  sway  force  is  generated 
for  this  configuration. 

3.   Rudder  and  Propulsion  Sway  Forces 

They  are  determined  from  equation  (11-11  (b) )  and 
(11-12) . 

The  sway  forces  due  to  the  seals ,  air  plenum  and 
aerodynamic  drag  are  neglected  in  this  modeling. 
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I.   YAW  MOMENTS 

1.   Sidewall  Yaw  Moment 

The  sidewall  yaw  dynamic  equations  are  given  in  [4] . 
After  some  simplification,  the  net  widewall  yaw  moment  is 
calculated  from 


N=«  =  N,  +  N_  +  (X^w  -  X   )  b         (11-44  (a)) 
sw      1    2     pw    sw   o 


where 


U1     =   2u(v  +xsr  -  zsp)A22£  (1  -  X) 


-  7r/5-pUd  ±  bQ$)2lu(v  -  zsp)        (11-44  (b)) 


N2   =   "  P  Cd  id(2  s1  i3  r  V  +  s2  I2   Vr2)         (11-44  (c) ) 

I   +  X 
A  =  — *■ — -  =  The  ratio  of  C.G.  location  of    (11-44  (d) ) 

sidewall  bow  to  the  craft 
length 

s,  =   1-3M1-X)  (II-44(e)) 

s2   =  2X2~   1  (11-44 (f)) 

Xpw  =   ~  V2  p  C  Aw  u2  (11-44  (g)) 

*  P 

=  port  wetted  area  surge  resistive  force 


Xsw  =   -  1/2  p  Cf  Aw  u2  (11-44  (h)) 

=   starbord  wetted  area  surge  resistive 
force 
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Aw  =   [2Ud  ±    bo*]    +   bsH  (11-45) 


fa 

s 


6£  +  bd 


=  Average  sidewall  flat  surface  width  at  the  keel. 

In  equation  (11-45)  ±  term  related  to  the  starboard 
or  port  wetted  surface  calculations. 

2.   Rudder  and  Propulsion  Yaw  Moments 

The  yaw  moments  due  to  the  rudder  and  propulsion  sys- 
tem are  determined  from  equation  (11-11  (d) )  and  (11-13 (b) ) 
respectively. 

The  seals,  air  plenum  and  aerodynamic  yaw  moments  are 
neglected.   Actually  the  aerodynamic  yaw  moment  should  be  taken 
into  account.   Since  during  the  turn  maneuver,  the  yaw  drift 
angle  is  produced,  thus  a  positive  aerodynamic  yaw  moment 
contribution  to  the  yaw  motion  has  been  generated.   From  6 
DOF  nonlinear  model  simulations ,  this  contribution  is  determined 
at  5.2%  of  the  total  yaw  moment. 

J.   TOTAL  PITCH  MOMENT 

The  summation  of  the  all  dynamic  pitch  moments  gives  the 
net  pitch  moment  which  can  be  written  as 

M   =   M  c    +  M    +  M.  -  +  M.   +M..+M    +M       (11-46) 
sf    sa    of    Da    pi    sw    p 
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K.   TOTAL  ROLL  MOMENT 

Before  calculating  net  roll  moment,  the  roll  moment 
equation  due  to  the  cross  flow  drag  term  must  be  derived. 
Including  also  the  added  mass  effect  to  the  equation  (11-40) 
and  (11-41) ,  the  result  is 


K,   =   (D   +  Y    ) (b   cos  6  -  z   sin  6)     (11-47) 
d       sw    sw.    o  o 


Combining  all  roll  dynamic  equations,  the  net  roll 
moment  is  calculated  as 


K      =      K    -    +    K        +   K,  .    +   K,        +    K         +K         +K+K+K, 
sf  sa  bf  ba  sw  pr  r  p  d 


L.   TOTAL  SWAY  FORCE 

From  equation  (Il-ll(a)),  (11-12),  (11-42)  and  (11-43), 
the  total  sway  force  is  obtained  as 


Y=Y+Y    +Y    +Y  (11-48) 

r    pr    sw    c 


M.   TOTAL  YAW  MOMENT 

The  summation  of  the  sidewall,  rudder  and  propulsion 
yaw  moments  gives  the  net  turning  moment  about  the  z-axis 


N   =   N    +  N   +  N  (11-49) 

sw    r    pr 
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N.   FIVE  DEGREE-OF-FREEDOM  EQUATIONS  OF  MOTION 

Considering  equations  of  motion  for  XR-3  test  craft, 
the  yaw  (i|i)  is  defined  as  the  rotational  motion  of  the  craft 
coordinate  system  relative  to  the  fixed  frame  and  the  positive 
if)  angle  is  measured  in  the  clockwise  direction.   This  sign 
convention  is  in  agreement  with  the  rudder  order  sign  defini- 
tion.  Since  positive  rudder  angle  corresponds  to  turning 
to  the  left,  producing  negative  yaw  moment  and  finally  at 
the  equilibrium  condition,  the  negative  yaw  moment  is 
generated. 

Assuming  small  values  of  pitch  (8)  and  roll  angles  and 
relating  heave  motion  to  the  free  surface,  the  relationship 
between  the  fixed  and  craft  coordinate  system  can  be  written. 
Eliminating  the  terms  related  to  the  fixed  coordinate  system, 
the  simplified  5  DOF  equations  of  motion  are  obtained  as 
follows; 

m  v  =   Y  (11-50) 

m  w  =   Z  (11-51) 


I   p  -  I    r  =   K  (11-52) 

xx  c  xz 


I    q   =   M  (11-53) 

yy 


I    r  -  I   p  =   N  (11-54) 

zz      xz  c 
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The  right  hand  sides  of  the  equations  of  motion  give  the 
forces  and  moments  acting  on  the  craft  in  terms  of  the  craft 
in  terms  of  the  craft  coordinate  system  variables.   In 
equation  (11-57)  w  represents  heave  acceleration  and  all 
other  variables  in  these  equations  have  already  been  defined. 

In  order  to  develop  a  5  DOF  linear  model,  only  the  surge 
equation  of  motion  (x-direction)  is  neglected.   Since  from 
6  DOF  nonlinear  simulation  results,  very  small  velocity 
reduction  is  observed  during  the  turn  maneuver  of  the  XR-3 
test  craft  (over  the  speed  range  20  through  30  knots) . 
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III.   LINEARIZATION  OF  THE  EQUATIONS 

A.   TAYLOR  SERIES  EXPANSION  AND  TOTAL  DIFFERENTIAL 

Generally,  the  linearization  process  can  be  performed  by 
means  of  a  Taylor  series  expansion  using  the  linear  part  of 
the  series  (omitting  high  order  terms) .   The  same  process  can 
be  easily  done  applying  a  total  differential  to  the  net  force 
and  moment  equations,  since  the  effect  of  small  disturbances 
to  the  system  is  to  be  determined. 

The  total  differential  of  a  given  dependent  variable  A 
with  respect  to  the  independent  variables  x,  y  and  z  is  given 
as 


j,      3A  ,   ,  3A  ,   .  3A  ,  ittt    i \ 

dA  =  -r—  dx  +  tt—  dy   +  it—  dz  (III-l) 

3x      9y  -*    dz 


B.   LINEARIZATION  OF  THE  EQUATIONS  OF  MOTION  AND  AIR  FLOW 
DYNAMICS 

The  total  differentials  of  the  equations  of  motion  and 

air  flow  dynamics  have  been  derived  as 

dw  =  dZ/m  CIH-2) 


dq  =   dM/I  (III-3) 

yy 


Pa    1/2   _ 

dm.   =   -[pa  n  kD  +  C  A.C — )X/ZJdP,     (III-4) 

*  a    p    n  I   2-  (Q)       b 

a 
dv   =   dY/m  (III-5) 
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where  Ph(0)  represents  the  steady  state  value  of  the  plenum 
gauge  pressure  according  to  the  final  equilibrium  condition. 

The  equations  (11-52)  and  (11-54)  both  involved  second 
order  yaw  and  roll  acceleration  terms.   Therefore  applying 
the  total  differential  and  solving  these  two  equations 
for  dp  and  dr  we  have 


dr  =   dK/I  +  dN/I  (III-6) 


dp  =   dK/I,  +  dN/I  (III-7) 


I,  I  and  I,  are  defined  as 
r      k 


A   *     *     "  I   2 

I  £        xx   2Z **-  (III-8(a)) 

xz 


I   4*   xx  zz ^-  (III-8(b)) 

r         I 

xx 


A    I     *     "  I    2 

Ik  *»  -^ — 12 xz_  (IH-8(c)) 

zz 


Equation  (III-6)  and  equation  (III-7)  contain  both  net 
incremental  roll  (dK)  and  yaw  (dN)  moments.   This  is  simply 
a  cross  product  inertial  term  (I   )  effect  in  the  equations 
of  motion. 

In  equation  (11-59),  dP,  is  not  a  state  variable,  it  can 
be  considered  as  an  output  variable,  therefore  from  equation 
(11-26)  and  (11-27),  dP,  must  be 
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3Pb       3Pb 

dpb  =  -W  dz  +  aiT  dmb 


C.   FIVE  DEGREE-OF-FREEDOM  LINEAR  MODEL  SENSITIVITY  COEFFICIENTS 

The  sensitivity  coefficients  related  to  the  right  hand  side 
of  the  equations  of  motion  are  defined  as  the  partial  deriva- 
tives of  the  net  forces  and  moments  divided  by  craft  mass  or 
inertial  terms. 

In  [2],  Barnes  has  developed  a  2  DOF  (heave  and  pitch) 
linear  XR-3  model.   Therefore  most  of  the  neave  and  pitch 
sensitivity  coefficients  are  directly  used  in  order  to  develop 
a  5  DOF  linear  model.   However  some  modification  is  made  for 
plenum  sensitivity  coefficients  and  pitch  mode.   Also  some 
of  the  steady  state  operating  point  values  have  been  reinves- 
tigated using  the  6  DOF  nonlinear  model. 

In  all  the  linearized  equations,  the  final  steady  state 
values  of  the  variables  have  been  used. 

1.   Pitch  and  Heave  Sensitivity  Coefficients 

In  [2]  equation  (111-40) ,  the  sensitivity  coefficient 
due  to  z  has  been  redefined  as 


9P,    _     fa   £q 
DHPZ   =   -  Ab(0)  j^  +   Pb(0)  -Sjj-2- 


This  is  claculated  from  (11-22)  where  the  second  term 
does  not  appear  in  [2] . 

The  modified  form  of  the  plenum  pitch  sensitivity 
coefficient  due  to  state  variable  z  is  to  be 
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DPPBZ   =   -  Zp(0)  ^-Xcp(0)  Pb(0)  JSp. 


The  derivation  of  these  modified  equations  is  based 
upon  the  total  differential  of  the  plenum  lift  force,  it 
is 


3Z_        3Z 
dZ 


P 


i£d*»  +  <ih» 


From  equations  (II-9) ,  (11-26)  and  (11-27),  the  total 
differentials  dV,  and  dP,  have  been  calculated  as 

fa  I 
dVb  =   ~(Ab<°>  "  Zd{0)      \   }  dz      CIII-9) 


dpb  '  v   pb(0)[M^ToT  +  v^ToT(Ab(0)-^d(0)  -V^)]d2-   (III"10) 

dP.  and  dV,  equations  are  also  different  from  the  equations 
derived  in  [2] . 

The  bow  and  aft  seal  pitch  moments  sensitivity 
coefficients  (DPSFTH  and  DPSATH)  due  to  6  and  plenum  air 
mass  equation  (AMB)  are  rederived  and  their  modified  forms 
are  used  for  the  pitch  mode  in  the  5  DOF  linear  model  computer 
program. 

Considering  the  5  DOF  linear  model  state  space  matrix 
representation,  the  a..'s  represent  "A"  matrix  elements. 
Therefore  the  a. .  values  relating  to  the  sensitivity  coeffi- 
cients for  heave  motion  are  determined  as 
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a21   =   DZZ 


a23   =   DZTH 


a0yl   =   DZDP 
24 


a25   =   DZMB 


where  DZZ,  DZTH,  etc.,  have  been  derived  in  [2].   However 
their  modified  forms  can  be  found  in  the  5  DOF  linear  computer 
program. 

For  pitch  motion,  a. .  coefficients  are 


a..   =   DTHZ 
41 


a43   =   DTHTH 


a..   =   DTHDTH 
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a45   =   DTHMB 


2.   Plenum  Air  Flow  Sensitivity  Coefficients 

These  coefficients  are  calculated  directly  from  [2] . 
The  a. .'s  related  to  plenum  air  sensitivity  coefficients  are 


a5,   =   DMBZ 
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a55   =   DMBMB 


3.   Sway  Sensitivity  Coefficients 

Before  deriving  the  sensitivity  coefficients  for 
sway  motion,  some  partial  derivatives  related  to  the  rudder 
and  sidewall  hydrodynamics  must  be  calculated.   From  equation 
(11—11 (b) )  Y   is  already  known,  this  follows  the  partial 
derivatives  of  Y  with  respect  to  the  state  variables 

3Y 
Rp  *  3T  =  Ct(Ce/u)  Zs  (Ill-ID 

9Y 

R    =   *■—  =   -  C.  (C  /u)  (111-12) 

v     3v        t  e 

3Y 

R   =  -Kzr-     =      -  C.  (C  /u)  X  (111-13) 

r     5r         t  e     r 


From  equation  (11-39) ,  the  sidewall  cross  flow  drag 
terms  D   partials  are  derived  in  the  following  equations. 
Let  us  define  first  some  coefficients 


dl  "  1   P  Cd  l 


dr     "      Zd{0)    +  si9n<vr>    fa0   <t>(0) 


The  partials  are 
3D 


D 


sw       ,     .  ,„  x  TI  2 


z       3z 


d1  sign(Vr)  Vr"(0)         (111-14) 
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D.   =   -  d,  V  2(0)  b  (111-15) 

<p        1  r      o 


D   =   2  d,  d   sign(V  )  V  (0)(Z  -£,(0)/2)   (111-16) 
p        1  r       r   r     s   d 


D   =   -  2  dn  d  V  (0)  (111-17) 

v  1  r  r 


D   =   -  2  dn  d   sign(V  )  V  (0)  X  (111-18) 

r         1  r   *   r   r     s 


d2  and  d.,  are  defined  as 


d2   =   2  A33  bQ2   u  (111-19) 


d3   "  A22  U  (111-20) 


where  u  is  the  craft  speed  assumed  to  be  constant. 
For  sway  mode ,  the  state  matrix  elements  are 


a96   =   £<!?>   =   E%  Sin  6)  (III"21) 


where  D,  is  given  from  equation  (111-15) . 


1tdY\ 
a97   "  ml3p' 

=  |[R  +  (R  +d3Zg)sin  6  +  d3Zs]       (111-22) 
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a99     m*3v' 


i[Rv  +  (Dv-d3)sin  6  -  d3]      (111-23) 


a911     mlIrJ 

=  m[Rr  +  (Dr"d3xs)sin  6-  d3Xs-mu] 

(111-24) 
4.   Roll  and  Yaw  Sensitivity  Coefficients 

Recalling  equations  (III-6)  and  (III-7) ,  the  net 
roll  and  yaw  moments  appear  in  both  equations  of  motion.   In 

4 

this  case,  to  determine  sensitivity  coefficients  is  not  easy, 
since  the  yaw  and  roll  partials  are  involved  in  the  dynamic 
equation  with  different  inertial  terms  and  coefficients. 

In  order  to  simplify  the  problem,  the  roll  and  yaw 
partial  derivatives  are  calculated  individually,  then  combining 
these  coefficients  and  dividing  by  corresponding  inertial 
terms  give  the  system  sensitivity  coefficients  for  the  roll 
and  yaw  modes . 

a.   Roll  Mode 

The  net  roll  partial  derivatives  are  defined  by 
k.  .  variables. 


The  partial  derivative  with  respect  to  z  is 


k71  '  If  "  1¥(Ksf  +  Ksa  +  Kbf  +  Kba)     (III"25) 


In  equation  (111-26)  all  roll  components  have 
already  been  derived  in  the  preceding  sections . 
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From  equations  (11-30)  and  (11-31) ,  the  seal 
partials  are 


3K  .        b    o 

"5«    =   "  (2T3)   p  ^1/sin  31<5)  ^0) 


9K  b    , 

-g—  -   -  (273)   P  gd/sin  32°)  <J>C0) 


For  buoyancy  partials,  let  us  define  the  following 
coefficients; 


b1  =   2  p  g  £1  (111-26) 


b2   =   2  p  g  %2  (111-27) 


b3  A    (^_»,  +  _f  (11I.28) 


b4  =  taHX  +  6f  <IJI-29) 


b6   '  SET^  +  "a  (III-30> 


From  equations  (11-35)  and  (11-36) ,  the  buoyancy 
partials  are 


9K,  -       b  b 

-W-     =      -tE^(b5+b3).*(0) 


3K         b  b 

,ja   -     "   '  (b,+b,)(K0) 


3z        tan  62  v  6    4 
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Equation  (11-26)  follows; 


9K  3P, 

jf-     =  -<V0)  Pb(0)  +(Zs-£d(0)  ^£]  4(0) 


Similarly  for  sidewall  cross  flow  drag  partial  is 


■*-=■  =  D  (fa  cos  3  -  Z  sin  8) 

dZ  Z    O  O 


9  are 


where 


The  roll  partial  derivatives  with  respect  to 


3K 
K73      36 


-      3lT(Ksf   +  Ksa   +  Kbf   +  Kba>  (III"31) 


3K  fa         ^ 

-gfi     =      (273)      P  gU3/sin  31°)    <M0) 


3K  fa        - 

sa  ,  w   x  3 

36 


=     -    (373)      P   gU4/sin  32°) 


3K,  .  b     b 

"#     =      tlHH^S   +   b3    V    *(0) 


3K  b     b 

-TT    '     "  tlr|j(b6  +  b4  V    «>(0) 
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IS 


The  roll  partial  derivative  with  respect  to  M, 


3K 


'75 


)Mb 


3P. 


=  V0)(zs-£d(0))3Mr*(0) 

b 


(111-32) 


The   roll  partials  with  respect  to   <j>   are 


"76 


3K 

34> 


where 


3  4> 


(K    ,   +   K        +   K,  .   +    K,       +   K      +  K.) 
sf  sa  bf  t>a  p  d 


(111-33) 


3K 


sf 


34> 


-    (bw/2.3)J    p    g   Xi 


3K 


sa 


34> 


=      -    (.b^/2. 3)J    p    g 


3K. 


bf 


34> 


-  bx   b3      b5 


3K 


ba 


34> 


"   b2   b4      b6 


3K« 

jf     =     V0)     (Zs   "   V0))    l    bw 
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where 


where 


3K 

W~     ■  V6o  cos  e  "  Zo  sin  e) 


The  roll  partials  with  respect  to  p  are 


*77 


3K 

3P 


^2-(K    +  K,  +  K  ) 
3p   sw    d    r 


(111-34) 


3K 


sw 


1p 


-  d,  -  7.J   A„  u  -  32/tt  fa  2  B 
2     s    22  o    c 


3K 

c 

3p 


3Y 


sw. 


(Dp  +  ~lp^)     (fao  cos  6  -  Zo  sin  6) 


The  roll  partials  with  respect  to  V  are 


•79 


3V 


■577  (K    +  K ,  +  K  ) 
3  V    sw     d     r 


(111-35) 


9K. 
c 

3V 


(Dv-d3)  (bQ  cos  B 


Z   sin  6)  +  d,  X   Z 
o  3   s   s 


3K 
sw 

~5v~ 


d_  Z 
3   s 
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9K 

2 

w 


CMC  /u)  Z 
t  e     r 


The  roll  partials  with  respect  to  yaw  rate  are 


'711 


3K 

3r 


3r-(Ksw  +  Kd  +  Kr) 


(111-36) 


where 


3K 


sw 


3r 


d-  X   Z 
3   s   s 


3K 


sw 


9r 


(Dr  -  d3Xg) (bQ   cos  6  -  ZQ  sin  B) 


3K 

ir-Z-     =   C.  (C  /u)  X   z 
3r       t  e     r   r 


b.   Yaw  Mode 

The  following  coefficients  are  defined  for  yaw 
rotational  motion; 


nl       A22  £(1  "  A) 


(111-37) 


n2   =    tt/5  p  u  £ 


(111-38) 


n3  =  -  sign(Vr)  p  Cd  &d(0) 


(111-39) 
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Cx   =   2  S1    I3 


(111-40) 


c2  A  s2  s2 


(111-41) 


Xw   =   -  1/2  p  Cf  u2  I 


(111-42) 


The  net  yaw  partial  derivatives  are  defined  by 


n. ■  variables. 


The  yaw  partials  with  respect  to   <J>   are 


3N 


n 


sw 


116 


8<j> 
3N 


cT$ 


d(j>       pw 


x     )    b 
sw       o 


(111-43) 


where 


3N 
9^" 


4   n2   V(0)    bQ*    (J>(0) 


^2-(X        -    X      )       =      -    4    X 
9<j>      pw  sw  w 


The  yaw  partials  with  respect  to  p  are 


n 


117 


3p  dp     sw  r 


(111-44) 


The  N   and  N  partial  derivatives  can  be  calculated  from 
equations  (11-44  (a) )  and  (11-11 (d) ) . 
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The  yaw  moments  partial  derivatives  with  respect 


to  V  and  r  are  derived  similarly,  the  n. .  coefficients 
related  to  these  modes  are 


n 


119 


9N 
9V 


(111-45) 


n 


9N 
1111  "      9r 


(111-46) 


c.   Roll  Sensitivity  Coefficients 

The  roll  sensitivity  coefficients  corresponding 
to  state  matrix  elements  a. .'s  are  determined  from  equation 
(III-7) .   They  are 


a71  "   k71/:Ck 


a73   "   k73/Xk 


a75   =   k767lk 


a76   =  k76/:Ck  +  n116/:n 


a77  =  k797lk  +  n119/:n 


a711  =  k711/:Ek  +  nllll/:E 


where  inertial  terms  I  and  I,  are  found  from  equations  (III- 
8(a))  and  (III-8 (c) ) . 
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d.   Yaw  Sensitivity  Coefficients 

For  yaw  motion,  the  sensitivity  coefficients 
a. .'s  have  been  determined  from  equation  (III-6) . 
These  coefficients  are  listed  below; 


alll  "  k71/:L 


a113   "   k73/3: 


a115   "   k75/:E 


a116   =   k76/:E  +  n116/:Er 


a117   =   k79/:C  +  n119/:Lr 


allll   =   k711/:E  +  "llll^r 


where  the  inertial  term  I   is  calculated  from  equation 
(III-8(b)). 

D.   STATE  SPACE  REPRESENTATION  AND  OUTPUT  EQUATION 

The  state  vector  of  the  linear  5  DOF  model  is  shown  as 


X      = 


=      [dZ   dW  d9   dq  dJV^  d<f>   dp  dy   dV  dip   dr] 


(.111-47) 


E"hen,  the  system  state  equation  must  be 


X  =  Ax  +  Bu 


(111-48) 
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And  the  output  equation  is 


y  =  C  X  (111-49) 

"A"  matrix  elements  have  already  been  determined.   However 
for  "B"  matrix  elements,  the  step  input  to  the  system  can 
be  considered  as  a  weight  removal  disturbance  and  for  the 
roll  and  yaw  modes,  a  rudder  order  disturbance. 

For  weight  removal  forcing  function,  the  input  is  to  be 


W  -  W. 
b21   =   -V  (III-50) 


where 


Wf  =   final  craft  weight 

W.   =  initial  craft  weight 

m   =  craft  mass  corresponding  to  the  final 
weight  W_. 


When  the  rudder  order  is  used,  sway,  roll  and  yaw  forcing 
functions  are  generated  as  follows; 


sway        b91  =  i(2  Xt  sin  6r  +  Ct  V         t111"51* 


roll        b.,,   =   (2  X,_  sin  5   +  C.  6  )  (-  z  /I,  +  x  /I) 
71         t      r    t  r     r  k    r 


yaw         b, ,,   =   (2  X,  sin  6  + c.  5  ) (x  /I  +  z  /I) 
j  111         t      r   t  r   r  r    r 
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Finally,  defining  state  variables  x,,  x_ ,  ...,  etc., 
the  state  equations  of  the  5  DOF  linear  model  are  given 
as  follows; 


Xl      X2 


X2   "   a21Xl  +  a23X3  +  a24X4  +  a25X5 


X3   ~   X4 


X4   ~   a41Xl  +  a43X3  +  a44X4  +  a45X5 


X5   "   a51Xl  +  a55X5 


X6  "      X7 


X7   "   a71Xl  +  a73X3  +  a75X5  +  a76X6  +  a77X7 


+  a7gx9  +  a711xn 


Xg      x9 


x9   "   a96X6  +  a97X7  +  a99X9  +  a911Xll 


x10    Xll 
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Xll   "   alllXl  +  a113X3  +  a115X5  +  a116X6 


+  a117x7  +  a119xg  +  a1111x11 


All  coefficients  of  the  state  equations  are  defined 
in  the  preceding  sections. 

Considering  outputs  of  the  linear  model,  only  incremental 
changes  of  the  variables  can  be  observed. 
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IV.   LINEAR  FIVE  DEGREE-OF-FREEDOM 
COMPUTER  PROGRAM  DESCRIPTIONS 

The  Digital  Simulation  Language  (DSL/360)  is  used  for 
the  linear  5  DOF  XR-3  computer  program.   Eleven  first  order 
integrator  blocks  are  needed  for  five  degree  of  freedom  simu- 
lations and  for  all  linear  runs,  fourth  order  Runge-Kutta 
with  fixed  interval  integration  method  is  selected. 

In  addition  to  the  craft  dimensional  values,  initial  and 
final  steady  state  operating  points  are  entered  into  the 
program  for  linearization  purposes  and  plot  outputs. 

The  input  parameters  of  this  program  are 

V    =  Craft  speed  (ft/sec) 

WF   =  Final  craft  weight  (lb) 

W    =  Initial  craft  weight  (lb) 

RUDANG  =   Rudder  input  (positive  for  port  turn) 

DR   =  Average  deadrise  angle 

PBBAR  =   Final  steady  State  value  of  plenum 
pressure  (lb/ft  ) 


ALD  =  Final  steady  state  craft  draft  (inches) 
PHI   =  Final  steady  state  roll  angle 


FYP  =  The  propulsion  system  thrust  corresponding 
to  the  craft  speed  (for  each  of  propellers) 

VL   =  Steady  state  lateral  speed  (ft/sec) 

R    =  Steady  state  yaw  rate  (rad/sec) . 
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The  steady  state  values  of  the  variables  are  determined 
from  6  DOF  nonlinear  model  simulations. 
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V.   SIX  DEGREE-OF-FREEDOM  NONLINEAR  MODEL  SIMULATIONS 


In  this  chapter  6  DOF  nonlinear  model  simulations  are 
discussed. 

At  first, some  of  the  runs  are  made  in  order  to  determine 
craft  initial  steady  state  operating  point  corresponding  to 
craft  speed  and  thrust  inputs.   Since  6  DOF  model  weight 
transient  response  can  be  analyzed  by  means  of  entering  proper 
craft  initial  values  into  the  program,  this  is  simply  an 
initial  condition  response  of  the  nonlinear  dynamic  system. 

A.   WEIGHT  REMOVAL  TRANSIENT 

For  30  and  20  knot  craft  speeds,  the  steady  state  values 
of  some  variables  are  listed  in  Table  I. 

Table  I 
Steady  State  Values  of  the  Weight  Removal  Response 


Speed 

Pitch 
(degree) 

Pb 

Draft 
(inch) 

Thrust 
(lb) 

Weight 

(knt) 

(lb/ft^) 

(lb) 

30 

0.95 

24.8 

8.6 

355.6 

6722 

30 

0.26 

24.8 

5.6 

287.8 

6050 

20 

0.45 

24.8 

6.06 

186.2 

6050 

For  30  knot  craft  speed,  two  different  steady  state  pitch 
angles  are  observed  corresponding  to  the  different  weights. 
In  [2]  for  pitch  transient  response,  0.65°  steady  state  pitch 
angle  value  has  been  used  for  both  linear  and  nonlinear 
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simulations.   On  the  other  hand  as  seen  in  Table  I,  0.26° 
and  0.95°  steady  state  pitch  values  are  obtained.   This  dis- 
agreement may  be  explained  by  noting  that  different  dimensional 
values  were  used  or  more  probably  a  different  modified  6 
DOF  computer  program  deck  was  used. 

For  6050  lb  craft  weight,  the  steady  state  pitch  angle 
0.26°  is  in  agreement  with  the  value  determined  by  Menzel  in 
figure  6  in  [6] . 

During  the  weight  removal  response,  simulations  of  6  DOF, 
two  different  stern  seal  base  leakage  areas  are  used.   In  [1] 
a  smaller  leakage  area  is  assumed  for  the  linear  heave  only 
model  as  A„  =  0.384.   However  Ag    =  3.84  also  has  been  simulated 
and  it  is  understood  that  the  smaller  leakage  area  is  more 
reasonable  considering  other  steady  state  variables  in  the 
air  flow  dynamic  equations  (11-23)  and  (11-24) .   Since  at 
the  steady  state  condition  Q.   =  Q   .  must  be  satisfied,  that 
permits  only  smaller  A,  for  a  given  Q  and  P,  steady  state 
values.   Therefore  A„  =  0.384  is  preferred  and  used  for 
nonlinear  simulations. 

B.   TURNING  MANEUVER 

When  the  3  DOF  pitch,  heave  and  roll  linear  model  was 
developed,  the  6  DOF  nonlinear  model  simulations  were  used 
in  order  to  determine  all  steady  state  roll  moments  during 
the  turn  maneuver.   In  addition  to  the  roll  moments,  yaw 
and  sway  outputs  were  also  investigated  and  the  results  are 
tabulated  in  Tables  II,  III  and  IV. 
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Table  II 
Steady  State  Roll  Moments 
(rudder  =  5°,  30  knot  speed) 

Roll  Moments Positive  Negative 

Stern  seal  FKSS  -  -2.2 

Bow  seal  FKBS  -  -28  0. 

Sidewall  FKSW  -  -41.5 

Rudder  FKRUD  131. 

Propulsion  FKP  0.057 

Aerodynamic   FKAED  -  -22.7 

Plenum  ABPB*PHI* (-Z)  215.3  - 


Table  III 
Steady  State  Yaw  Moments 
(rudder  =5°,  30  knot  speed) 


Yaw  Moments 


Stern  seal   FNBS 

Bow  seal   FNSS 

Sidewall   FNSW 

Rudder   FNRUD 

Propulsion  FNP 

Aerodynamic  FNAED  -  -  5  6.41 


Positive 

Negative 

3.22 

- 

0. 

0. 

- 

-482.2 

532.65 

- 

2.78 

— 
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Table  IV 

Steady  State  Sway  Forces 

(rudder  =  5°,  30  knot  speed) 

Sway  Forces  Positive  Negative 

Sidewall  FYSW  -              -351.8 

Rudder  FYRUD  -              -  47.6 

Aerodynamic   FYAED  -              -  19.07 

Centrifugal   -R*U*AM  419.24 

From  Table  II,  it  is  understood  that  for  a  given  (+) 
rudder  input,  the  sign  convention  of  the  craft  dynamics  was 
in  agreement  with  turning  to  the  left.   However,  a  negative 
sway  displacement  plot  is  observed  for  this  maneuver.   As  a 
result,  in  the  R.H.S.  subroutine  of  the  6  DOF  program,  the 
sign  correction  should  be  done  for  y-displacement. 

The  steady  state  sway  forces  acting  on  the  XR-3  test 
craft  in  a  turn  to  the  left  are  shown  in  figure  6.   As  seen 
in  this  figure  and  from  Table  IV,  the  rudder  sway  force  is 
in  the  steady  state  condition.   For  a  left  turn,  initially 
the  rudder  sway  force  is  generated  in  the  reverse  direction 
and  finally  a  negative  steady  state  rudder  sway  force  is 
produced.   Further  investigation  is  done  in  the  RUDDER  sub- 
routine dynamic  equations  of  the  6  DOF  computer  model ,  but 
no  error  has  been  determined  in  these  equations.   Also, 
different  craft  speeds  and  different  effective  rudder  angle 
coefficients  (ENDFAC)  are  used  but  every time  the  same 
configuration  is  observed. 
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Figure  6.   Steady  state  sway  forces  configuration 
at  the  left  turn 
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In  addition  to  rudder  dynamics,  the  propulsion  system 
of  the  XR-3  contributes  to  the  roll,  yaw  and  sway  motions. 
However  this  contribution  is  not  a  dominant  factor  compared 
with  the  rudder  effect. 

For  the  XR-3  test  craft  propulsion  system  configuration, 
the  direction  of  the  thrust  vector  is  also  changed  as  well 
as  the  rudder  deflection.   As  seen  from  Tables  II,  III  and 
the  simulations  with  different  rudder  order,  it  is  understood 
that  the  propulsion  thrust  component  due  to  the  rudder  order 
is  not  calculated  properly,  therefore  the  following  modifica- 
tion is  made  in  the  PROPULSION  subroutine  of  the  6  DOF  program, 

FYS   =   -STHSTS*CD+THSS*CD 

where  FYS  represents  the  sway  force  due  to  the  starboard 
propeller,  and  the  others  are  defined  in  the  propulsion 
subroutine  program. 

C.   DEADRISE  ANGLE  EFFECT 

As  reported  in  [4] ,  the  deadrise  angle  has  a  major  effect 
on  the  surface  effect  ship's  roll  behavior  and  therefore  is 
an  important  design  consideration.   The  deadrise  angle  effect 
on  the  roll  motion  has  been  investigated  by  Reidel  in  {3] . 
In  this  investigation,  some  arbitrary  coefficients  (PROMO  1, 
2,  etc.)  are  used  for  the  vertical  force  component  of  the 
sidewall  hydrodynamics.   From  this  investigation,  the  conclusion 
reported  by  Reidel  was  that  more  damping  could  be  obtained  for 
the  pitch  and  roll  motions  by  changing  the  deadrise  angle. 
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The  deadrise  effect  is  reinvestigated  by  simulating  the 
system  with  different  rudder  orders  including  zero  disturbance, 
then  the  heave  forces  due  to  sidewall  hydrodynamics  were 
compared  to  each  others*  As  a  result,  no  significant  sidewall 
heave  force  change  was  observed.   For  instance,  for  0°  and 
5°  rudder  inputs,  the  sidewall  heave  forces  (FZSW)  are  found 
to  be  -568.8  and  -568.05  respectively.   This  result  is  not 
acceptable,  since  in  a  turn  maneuver,  the  cross  flow  drag 
is  generated  at  the  outboard  sidewall,  so  the  deadrise  projec- 
tion of  cross  flow  drag  must  be  contributed  to  the  sidewall 
heave  force.   More  investigation  has  been  made  in  the  SIDEWALL 
subroutine  of  the  6  DOF  program  and  it  has  been  determined 
that  only  the  force  component  due  to  added  mass  given  by 
equation  (11-21)  is  projected  as  a  lift  force.   Actually  the 
cross  flow  drag  component  must  be  included  to  this  projection, 
since  the  frontal  area  for  cross  flow  drag  component  is  out- 
board of  the  submerged  section  of  the  sidewall.   Therefore 
it  is  concluded  that  the  XR-3  test  craft  6  DOF  program  SIDEWALL 
subroutine  does  not  include  effectively  the  lift  force  due 
to  deadrise.   On  the  other  hand,  the  lift  force  and  horizontal 
drag  component  can  have  different  drag  coefficients,  therefore 
the  modification  of  the  sidewall  subroutine  is  left  for  future 
studies. 
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VI.   COMPARISON  OF  NONLINEAR  AND  LINEAR  MODEL  RESULTS 

A.  C.G.  ACCELERATION  (see  figs.  9  and  23) 

A  10%  weight  removal  disturbance  is  used  for  both  models. 
The  linear  and  nonlinear  model  showed  very  rapid  response 
with  small  rise  time.   As  seen  in  figure  9,  the  6  DOF  model 
transient  is  relatively  more  damped  compared  with  the  linear 
model  response.   The  linear  model  showed  a  small  overshoot 
but  it  is  also  sufficiently  damped.   Small  oscillatory  behavior 
is  observed  in  both  figures  which  is  due  to  pitch-heave 
coupling.   In  general  the  two  models  agree  in  transient  response 
and  steady  state  behavior. 

Both  of  the  C.G.  acceleration  curves  were  normalized  by  a 
factor  "g'^the  gravitational  constant. 

B.  PITCH  TRANSIENT  (see  figs.  10  and  24) 

As  seen  in  these  figures,  both  models  showed  an  oscillatory 
transient  response  with  the  same  frequency,  w  =  3.12  rad/sec. 
However  their  steady  state  values  are  different.   In  fact, 
the  pitch  moment  due  to  plenum  pressure  and  the  effective 
plenum  pressure  center  play  a  dominant  role  for  pitch  motion. 
In  [2]  Barnes  developed  plenum  modeling  including  effective 
plenum  roof  area  calculations  and  plenum  pressure  centroid 
shift  due  to  draft  variation.   Therefore  the  linear  model 
plenum  modeling  is  different  from  the  nonlinear  modeling, 
so  that  a  greater  plenum  moment  arm  is  calculated  in  the  linear 
model. 
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The  steady  state  pitch  angle  values  are  0.5°  and  0.75° 
for  nonlinear  and  linear  models  correspondingly. 

C.  PLENUM  PRESSURE  TRANSIENT  (see  figs.  7  and  21) 

As  seen  in  these  figures,  both  models  are  exactly  in 

agreement  with  each  other.   In  each  curve,  the  plenum  pressure 

2 
suddenly  drops  to  21.6  (lb/ft  )  then  increases  up  to  its 

steady  state  value  with  sufficient  damping.   This  is  expected 

behavior  because  the  plenum  volume  increases  when  10%  weight 

is  removed,  obviously  this  disturbance  causes  a  sudden  drop  in 

plenum  pressure.   As  seen  before,  pitch-plenum  pressure  coupling 

had  been  observed  during  the  exponential  increasing  phase  of 

the  plenum  pressure.   At  the  fifth  second  of  simulation  time, 

nonlinear  model  plenum  pressure  is  about  24.29  while  the 

linear  model  pressure  is  24.25. 

D.  DRAFT  TRANSIENT  (see  figs.  8  and  22) 

Excellent  agreement  was  achieved  between  linear  and  non- 
linear model  draft  responses.   As  expected,  the  draft  decreases 
exponentially  when  the  weight  is  removed  from  the  craft. 
The  draft  is  found  to  be  5.7  inches  for  the  linear  model, 
5.85  inches  for  the  nonlinear  model  at  the  fifth  second  of 
simulation  time. 

As  a  result,  both  models  exhibited  the  same  draft 
responses. 

E.  ROLL  RESPONSE 

As  mentioned  before,  at  first  a  3  DOF  linear  model  (pitch- 
heave-roll)  had  been  developed.   3  DOF  linear  and  6  DOF 
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nonlinear  models  roll  behavior  comparison  is  done  by  means 
of  a  step  input  roll  moment  disturbance  to  the  linear  model 
and  a  turning  command  to  the  nonlinear  model.   In  figure  13 
6  DOF  nonlinear  model  roll  angle  response  is  shown  for  a  5° 
rudder  order  (left  turn),  in  this  simulation  1.18°  outboard 
steady  state  roll  angle  is  generated.   From  the  R.H.S.  sub- 
routine of  the  6  DOF  program,  all  the  steady  state  roll  moment 
components  are  determined  as  seen  in  Table  II  and  the  sum  of 
these  disturbance  moments  applied  to  the  linear  model  roll 
mode  as  a  step  input.   Roll  behavior  is  analyzed  from  figure 
27.   As  expected,  the  linear  model  showed  a  greater  overshoot 
and  more  oscillatory  response  compared  with  the  nonlinear 
roll  transient.   Since  the  linear  model  has  been  disturbed 
by  a  step  roll  moment  corresponding  to  the  nonlinear  model 
steady  state  roll  moments.   On  the  other  hand,  its  steady 
state  roll  angle  is  found  to  be  1.28°.   The  roll  steady  state 
difference  of  two  models  is  about  8.4%.   For  a  3°  rudder  input, 
nonlinear  and  linear  simulations  are  obtained  in  figures  11 
and  25  with  the  steady  state  values  0.62°  and  0.65°  respect 
tively.   Hence  the  difference  is  4.8%  in  the  steady  state 
values  for  these  simulations. 

Figures  12  and  14  show  XR-3  yaw  rate  plots  and  corres- 
ponding linear  simulations  are  shown  in  figures  26  and  28. 

As  a  result,  the  3  DOF  linear  model  and  6  DOF  program 
roll  simulations  are  in  agreement  with  each  other  within 
10%  difference. 
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As  discussed  in  the  previous  chapter,  some  problems 
are  observed  in  the  6  DOF  nonlinear  program.   Therefore 
developing  two  more  modes  for  the  3  DOF  model,  the  5  DOF 
linear  model  has  been  generated.   By  means  of  this  model, 
roll  behavior  of  the  XR-3  in  a  turn,  deadrise  angle  and 
cross  product  inertial  term  effects  have  been  investigated 
for  20  knot  craft  speed  with  different  rudder  deflections. 

The  roll  response  can  be  seen  in  figure  29  omitting  the 
cross  product  inertial  term  (Iv_  =  -2800.)  in  the  equations 
of  motion,  applying  a  10°  rudder  order  and  assuming  70  degree 
average  deadrise  angle  of  the  craft.   For  this  particular 
test,  the  craft  starts  inboard  rolling  with  the  oscillatory 
transient  then  goes  to  0.47°  steady  state  roll.   Using  the 
same  input  parameters  and  also  including  the  cross  product 
term  I   in  the  equations  of  motion,  smaller  overshoot  (18%) 
and  a  more  damped  roll  response  has  been  observed  in  figure 
31.   As  seen  in  this  figure,  at  the  initial  phase,  secondary 
overshoot  has  been  generated.   The  steady  state  roll  angle  is 
0.40°  and  the  oscillation  frequency  w  =  3.5  rad/sec. ,  for 
this  simulation.   On  the  other  hand,  the  corresponding  6  DOF 
nonlinear  simulation  is  shown  in  figure  17  and  also  a  roll 
rate  plot  is  shown  in  figure  18.   The  nonlinear  simulation 
in  figure  17  shows  a  more  damped  transient  but  greater  steady 
state  roll  behavior  (1.32°)  compared  with  the  corresponding 
linear  simulation  shown  in  figure  31. 
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Figures  30  and  32  present  linear  simulation  curves  for 
5  and  15  degree  rudder  inputs  respectively.   Both  curves 
exhibit  the  same  linear  model  response  discussed  before. 
Finally  the  linear  and  nonlinear  model  steady  state  responses 
are  summarized  as  follows. 


6  DOF 

5  DOF 

Rudder 

Nonlinear 

Linear 

Reference 

5 

0.57 

0.24 

figures  15 

and  30 

10 

1.32 

0.40 

figures  17 

and  31 

15 

2.2 

0.45 

figures  19 

and  32 

The  deadrise  angle  effect  is  investigated  by  using  differ- 
ent deadrise  angle  values  for  10  degree  rudder  input  and  2  0 
knot  craft  speed  simulations.   The  results  are  tabulated  in 
Table  V. 


Table  V 
Sidewall  Deadrise  Angle  Effect 


(rudde 

.r  =  10°,  20 

knot  speed) 

Deadri 
angle 

se 
(degree) 

Roll  angle 
(degree) 

Lateral 

Speed 

(ft/sec) 

Yaw  rate 
(rad/sec) 

45 

-1.01 

4.18 

-0.064 

50 

-0.84 

4.04 

-0.069 

55 

-0.62 

3.90 

-0.076 

60 

-0.30 

3.83 

-0.080 

65 

0.02 

3.73 

-0.085 

70 

0.40 

3.60 

-0.100 
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For  the  linear  model,  the  lateral  velocity  and  yaw  rate 
plots  are  shown  in  figures  33  and  34.   Both  of  them  showed 
sufficiently  damped  transient  response.   From  figure  33, 
the  lateral  speed  steady  state  value  is  found  3.6  (ft/sec) 
and  as  expected  for  turning  to  the  left,  the  negative  yaw 
rate  is  generated  as  -0.1  ( r ad/sec ) . 
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Figure  7.   6  DOF  Plenum  pressure  transient 
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Figure  9.   6  DOF  C.G.  Acceleration 
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Figure  10.   6  DOF  pitch  angle  transient 
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Figure  11.   6  DOF  Roll  angle  plot,  6   =30 
30  knot  speed  r 
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Figure  12.   6  DOF  roll  rate  plot, 
30  knot  speed 
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Figure  13.   6  DOF  roll  angle  plot, 
30  knot  speed 
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6  DOF  roll  rate  plot, 
30  knot  speed 
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Figure   15. 


6  DOF  roll  angle  plot, 
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Figure  17.   6  DOF  roll  angle  plot, 
20  knot  speed 
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Figure  18. 
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Figure  19.   6  DOF  roll  angle  plot, 
20  knot  speed 
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Figure  20.   6  DOF  roll  rate  plot, 
20  knot  speed 


RD=  1  5 


6   =  15° 
r 


90 


o 
o 


PLOT    IS    PLENUM    PRESSURE 
XR-3    30    KNOTS    SPEED 


en 


o 
o 

CVJ 


o 
o 

o. 


o 
o 


o 
o 


o 
o 


in 


o 
o 


o, 


0.00 


1.00 


2.00 


3.00 


4.00 


5.00 


XSCRLE=  1.00 
YSCRLE=  5.00 


UNITS/INCH  RUN   NO.  1 
UNITS/INCH  PLOT  NO.  1 
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Figure  22.   Linear  5  DOF  draft  transient 
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Figure  23.   Linear  5  DOF  C.G.  acceleration 
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Figure  24.   Linear  5  DOF  pitch  angle  transient 
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Figure  26.   Linear  3  DOF  roll  rate  plot,  6   =  3( 
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Figure  27.   Linear  3  DOF  roll  angle  plot,  6^  =   5' 
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Figure  28.   Linear  3  DOF  roll  rate  plot, 
30  knot  speed 
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Figure  30.   Linear  5  DOF  roll  angle  plot,  6r  =  5' 
20  knot  speed 
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Figure  31.   Linear  5  DOF  roll  angle  plot,  6   =10' 
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Figure  33.   Linear  5  DOF  lateral  velocity  plot 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

In  the  preceding  sections,  the  5  DOF  linear  XR-3  model 
was  developed,  then  weight  removal  response  and  roll  behavior 
are  compared  with  6  DOF  nonlinear  simulations.   Also  some 
of  the  craft  sensitivity  parameters  are  investigated  for  the 
XR-3  test  craft. 

The  results  of  this  study  are  summarized  and  recommenda- 
tions are  made  as  follows. 

1.  The  system  response  is  investigated  by  both  a  linear 
model  and  a  6  DOF  nonlinear  model  by  means  of  sudden  weight 
removal  from  the  C.G.  and  by  sway  force  disturbance  in  a 
turn. 

Comparison  of  results  shows  that  the  difference  is 
within  10%  for  both  models.   In  fact,  assuming  accuracy  in 
measurements  of  about  10%,  the  linearized  model  gives 
acceptable  results  as  compared  with  the  nonlinear  model. 
Obviously  the  linear  model  simulations  take  less  computer 
CPU  time. 

2.  Concerning  roll  response  during  the  turn  maneuver, 
the  product  inertia  term  (I   )  affects  transient  response 

XZ 

such  as  overshoot,  damping,  etc.   More  overshoot  and 
oscillatory  response  have  been  observed  without  I 
However,  this  term  is  found  to  be  -2  800.  from  the  6  DOF 
program.   It  is  suggested  that  this  number  be  verified. 
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3.  The  importance  of  deadrise  is  observed  for  roll 
steady  state  behavior.   The  XR-3  test  craft  steady  state 
roll  angle  depends  upon  the  average  deadrise  angle  of  the 
sidewall.   Different  average  deadrise  angles  have  been 
simulated  by  means  of  the  linear  model  and  it  is  observed 
that  the  craft  rolls  inboard  or  outboard  according  to  the 
sidewall  deadrise  configuration. 

The  6  DOF  nonlinear  model  SIDEWALL  SUBROUTINE 
must  be  modified  by  including  a  cross  flow  drag  component 
deadrise  projection  (lift  force) .   However  a  new  drag 
coefficient  is  needed  for  this  lift  component. 

4.  The  aerodynamic  effect  of  the  craft  contributes 
about  5%  of  the  total  yaw  moment. 

5.  The  XR-3  failure  mode  operation  (such  as  fan  failure) 
can  be  simulated  by  using  either  the  linear  or  nonlinear 
model.   But,  the  nonlinear  model  is  preferable  for  failure 
mode  simulations,   since  it  includes  sidewall  and  stern 
seal  leakage  configurations. 
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APPENDIX  A.   COMPUTER  PROGRAM  LISTING 


OOOOOO 

OO 

<\jm«i-tr»vOf*- 

>TfM 

i-li-lf-lt-tr-lf-l 

>fr>r 

oooooo 

OO 

oooooo 

OO 

«/)t/)«/)U)  </)</> 

V)V> 

u/U'uuuju.uj 

IUIU 

(/IVH/.'l/)V)V) 

</)«/> 

ooooooo 

OOOOOOO 

ooooooc 

«/)VJ  «/)«/)  I/)  «/)«/> 

U-IU-'U.  U,  U.  U'LU 
V></>l/>«/>«/>t/>«/» 


OOOOOO 

ooo>Oi-HCjm 
c\jc\jrqmmr<> 
OOOOOO 
OOOOOO 
«/)«/>  1/5  V?  VI  lO 


ooo 

en  ton 
coo 
ooo 

U-'U    u. 

vtvun 


UJ 

o 

O 

2: 


or 

< 

UJ 


U- 

o 
o 
to 


f*1 
f 


to 


z 
o 
o 


I- II  o 

XOUh  hi 

u.  »_JO  </) 

woa — i  > 

^  II  ZQ.  I 

<£<    z-j    a. 

or    < 

O^U4t-«-< 
UIUUOI- 
l-OI-Zt-i 


«rom    o 

•   »i-l      o 

CM  II  i-li-l    »0 
.<t_l    «CM  II 

<Njs:r\jm  ii  < 
nis  it  ii  oo 

II«MII 

UOOlQuQCOX 


to 

z 
o 

00      M 

r-    «/) 


a: 


in 

•  ■*/ 
i-i>ru> 
II    •<-» 
mi— o  • 

•  O  II  Oui 
c\lIU»-«  II 
II  ma.  ii  2 

WIDOWS 
i"*ICQXX> 


OO  lO 

O  •    •  C\1>*J 

•  <NJc\4r-  (M    r-  • 

©rt|»0lM%0>%^«MrHOU"\rMU\         #v> 
i-H    IU>     •(Mf-HCNIM     •    •     •     •     •        lf\rt 

MO   »f*    »_J— I    •  P-r-ICOOOtPt-l  || 
»iH^H«D«THiO>0  II    II      •  II  _l 

ii  qcq.03>0|-  ii  ii  ii  it  ii  ii  ii  ii  ii  noooa3 
:*<£<t^«/>  ii  cMOr-KMrn^rirv^r-aoi-icMi-itM  n  a.a. 
</>_i_j_iOq;  ii  •— • t i i i i — i_j_ioca:«/)t/>z-J_J 


co  vo 

(M    •    »C0 

m*r  n  ii  • 

•  II  <U.fH    • 
IMXUJUJ  II  o 


*  #  * MMMOM#  *  * 


*** 


107 


ooooooo 

ooooooo 
ooooooo 

IUIMIWIU  UMU  111 


OOOO 

OOOO 
OOOO 

cocooooo 


oou 
eocvjm 
»3-coao 

HOO 

ooo 
00  OOOO 

U'U.  U-' 

oo«/;oo 


3 
* 

O 

X 
CO 

CO 

* 


CM 


oca 

J3I 

I  U.OI- 
-JU-l-iO 
3mJh 
U-Q*  -* 
-J -J-l* 
<l<l.*CC 

II  II  U.U. 
U.LL 

»-«•-•  II  H 

oosa. 

_j  CD  CD  CD 


a 


< 
+ 
cc 

0.0 


•    • 


00 

CO 


Ooroin 

l\J|—    •••!"» 

oioowtcp  • 
<  m  o  u*  tvj  >r  ■*■  c\j  o 
wm  •  ii  ii  won  ii 
<-»»  tf»xr-coc\ioo 
no  m  ?-x  ii  i  it  m 

>  cj  ai<t  <  •-• — >~i  < 


2 
O 
t-i 
oo 

Ul 

> 

z 
o 
u 


CO 
at 
o 

■a 
u. 


CO 

cc 

Ul 

> 
z 
o 
<-» 

u. 

—I 

O 


HO 
Q-vO 
•fcro 

OS. 

•MM 

NhOOO 
•^tfOCOCOC 

ooouu 

•    (11101111 
©CMCJCJO 

>o>~*  ♦* 
ro  ii  i-hnjcc.^ 

IIQaCOCUf- 
<J<l-»(_)  ||  Ul 
UI0C  II  II  <C0 
CJOi-i<m*-  II 
<ui  a:  educe 
aQUumu 


o 
a 

u> 

z 


< 

a. 
o 


< 
z 


a. 


UUsO 
•  cm 

»r  ("T        r-  • 

piNom        <tcj  »o 

r-  ii  ii   *<ir-\  •  •cvjua 

•  cc<ir\  •  ioono 

rn<i-  ii  mo  ii  ii  h-^o 

OICCUJO  II     I   MQO  II 
IICOX-J-I  II  XX  II  LL 


Q. 


cc 

LU 

u 
a. 
ui 

ex 

UJ 

a 
o 

CC 


II 
o 

z 
< 

cc 


»-* 

O-l 

=?> 

O  I 
UI(M 

a:* 
* 


ooc 

-100 
Ul  II 


a 

<r 
oo: 
<o 

(XU.I 

oo 

Ul* 

CO 

o*z 

•<< 

OCMt-U 

U->s.J.0C 

-*CJ»—  II 

II  -I  II  o 

oo<i«iz 
oo  ii  t-< 

<UUJQ 

«.i-c£ 


tn 


##* 


*#* 


#*#      ** 


108 


oo  oooooooo 

oooo  oo  u>  cr>  <r  i/*  cr*  cr  cr 

oo  oooooooo 

oo  oooooooo 

OOOO  OOOO  OOOO  OOOO  OOOO 

U.-UU  U.U1U  UU'U-'U'U-' 

«"O0  00</;OOOOV»000000 


oo 

o 

OO 

Os<j» 

rvj 

nJ-IO 

mo> 

o 

OO 

oo 

l-J 

r-li-< 

oo 

o 

OO 

oooo 

CO 

O000 

lilUJ 

u> 

U-'UJ 

oooo 

OO 

OOOO 

oo 

z 
o 

CO 

a. 


CvJ 
* 


O 

u. 

o 


< 
cc 

UJ 

Ql 

o 


u. 

z 


a 
< 
a: 

UJ 

a 

WW       •* 

qui  «m>o»-« 

•  #sr  «ooi 

II  II  II  II  II  II 

i_lf-4r-lf-IQ.I-< 

a.(-a._iu.uL 


a 
u 

o 

z 
< 

OO 

a. 

-i 

CD 

< 

oc 

> 

UJ 


o 

UJ 
UJ 


3 
<~> 
_l 
< 
o 


<-l(\J 

0C<X 


l-l- 

** 

« 

UIU:-* 
|-»-«-«IM 

_!_)<!< 

I  +  +  + 

uaoo 

_J_lr-U\J 
<l<O0</> 

II  II  Jt-S 
r-4(M  II  II 
CJOftOJ 
_J_IOOCO 


fHCVJ 
-J -I 
« 

GO 

au 

IX 
a.  a. 
** 

OO 

•   • 

CNJOJ 
II  II 
<-H«\J 

<■< 


X 
UQ. 
0.W 
OOO 

xo 
+  o 

0"v  — 

•—  < 

CM—  S 

UJ     — o  < 

C-hl  O 

OOOC0  >v 

-I  »x^>  o 

Q-CNJCDCO  • 

#H3\  •-< 

MNCDO  — 

Q-^O^-J  * 

*    »0«1  * 

OOUJUj<-hc\JC\J0>-J~  — • 

##XX-J-'#i-<<#  < 

t-<fvj|— I— <<*•-•*—  O. 

«W*  #ujuj>«—  u-oc 

-wm'fMlfl*   I  U-CO--  CO 

ZZJJXX<Oh<m  Q. 

hhi-h<<*  »uj  »a  i  xo    «- 

OOOO  I   +  IXQC<\J_J3a — I       * 

—uQhi-<v<ai*«<Q:ai 

NN-IJUQJJ>-<llfl*<t-> 

lx<i<'-<waj;  +  -uHiL# 
ft—"" "*•  3-*h-u-cj  I  <_>aucu< 

UUII    M    II    II  l_)_ICL  Ji*< 

i-n-««-*c\i<-icMX«jcaco  I  +X 
_£jC_j_j_j-jqc  ii  x<<z<icc 

||    II  <l<<<  |l  O  II  <*•  II  >Q.  II 
rO^UJUJUJUJOU-Q-  M  rH  It   II  00 
s^  ^  oooo  oooo_joo  coco  coca 
<<xx<j.<q_xx<i<j..>q.< 


UJ 

CO 


o 
o 

— * 

<(NI 

o* 

X* 

oco 


aroc 


« 


<\j<\i 


# 
— * 
cvj  co 

0£  00 
OOO 

—a 
z  + 

I-* 


oc<     «<■>.>», 

<♦  l-l-CMi-l 

CO—  UJLUOO 

CQCM  XX-J-JOO— *"* 

O I  hi-  «KCt-lC\J«lUJ 

**l  **  +  +  OOOOI—  0" 

0+  O0O0    •    »33U<0      > 

«~  .,-j  xxc\ic\j++x—  * 
ctrsi-j  +  +  v.*^— — t-  *  _jim 
<<— «l-l»-<  — ^jcvj*  •-<■*  * 
coi— - -l-xxxiocccooa.3* 
cocc*UJa-a-t-l-oox*oo3 

O.CTOX**  DO-'-'  +  nqw 

I  WII-3if't-<ZZD    •O?- 

0»  CC*  WOT3S«-JO#  * 

•-<*c*oo>->  +  •M-|~<t,*Ljoo 

(JO  »X  +  I  00>«,N.«-0Xn)> 
ww«ji.c3cjrvJf-<f-<rvj  +  xut.nnr 
•*•  nnu-i^ v)v>a aoa: *  «i/i 
Z  «£  •—  -J  «  3  js  _ i  _j  c\j — u->  #  r\j 
JJI-lQ-^  II  II  «--».<<oo  II    •  »c\j 

It  (J  II  II  f-HCM  ||  ||  ||  ||  _*</)OlM< 
Z  II  C\IO0O0O0C\Jr-li-IC\l  II  <M  II  II  II 
•-•  —1 0  OO  OO  O0  _l  _l  3  3  00  CM  t-i  CM  fO 

cj<i  co  oj  c_>  o  cu  co  co  aj  co  <i  o  o  (_> 


CM 

V. 
a 

-J 
< 
I 

O0 

rsl 

I 

a 
* 
oo 

X 

I 

_l 
II 


cc 


*** 


*♦* 


109 


OOCJO 

ooooo 

o 

oo 

OO 

OO 

sfir»>or^ 

(fi«j->or*-oo 

o> 

l-<CM 

m<a- 

m>o 

i-Hi-H  t-H«-4 

CMCNJCMCMCM 

CM 

nim 

men 

r-r- 

r_H-l.--i.-H 

P-lr-lr-lr-lF-l 

t-H 

f-it-i 

r-*r-4 

r-lr-l 

oooo 

ooooo 

o 

oo 

OO 

OO 

00000000 

OOOO  00  OOOO 

00 

OOOO 

OOOO 

vnn 

uju.iu.iu' 

UL-U    LL'U.IU. 

III 

U'l— ' 

IVIU 

1>'UJ 

U)V>V)Ui 

OOOO  OO  OOOO 

l/l 

oooo 

OOt/> 

«/>«/> 

oooo 

-hcmco-j- 

00  0000  00 


oooo 
00O-IO0O0 

U  lA'U-l." 
UlVHftVI 


a 
a. 

oc 

HI 

a. 
O 

-J 
<l 


oo  Hw 

_J  O-O-LU 

<  **oc 
3>  OOO. 
Q  #*X00 

»-.  oo  +  iu 

oo  XX<oC 

IJU  QCOCQ-CL 

r-4CM  +  + 

o  _j_ju-Z< 

2_  4      I-h           <<u.<00 

<  I-     X        gjujx_jac 

CM  UJ       CLt-n-iOOOO  +  O-  "♦■ 

►  J                                  X      *XXXX<.CLl_. 

OO  <I                                    »—      —  Q.U.  #*</)  +  OO 

|—  LU                                    *       D#*C)ClI<ia 

2  to  OO  5-  _lr-lCM*  #  + 1/>  + 
UJ  <  •  •  OC  <S3#*U.O.< 
Z  i-ICMi-l*                      CMfvJ      <        |  aJaj—*,»O0  +  CD 

O  O-OOOO-J— •                     ^V.      _J      00*  *  CMCMXU-C_00>-X 

£  tfO33<0                     i-tCM      Cv      MCMCMU"Mn  +  O0  +  O0>X 

•■  <x**ui  •«        _j_i    *     —  **   •  •<au.<»-H»-< 

CO  CO*i-*CMO0CM>-t-           <<U-0      ***CMCMC0  +  0_3_<« 

LU  COOOOO«+UJl_JLL      UJLUO0.J<O0r-l<M>»^X<QC<I-,v^ 

O  CLMJ_J_J*  OC-LXa3<l/>O0_IQ.00UJ_J  — 1X1  +  OJ  +  -^CM(T> 

ac  *  ot 4 4tt <lf-l—  a. cuj«;xu.ir  _j  acajou | — | —  u_ a. oor-i vi ooj- 

U  cull*  *«a.aj-t(-  *  *  J.  I   f  -H-    »IJ.CL**OOC_  +  UJ.OOU_UJ  | 

u.  <a.ir-i(\jcDuj(_)oiri*  oi-ru-iM*  x<-ii>J'-«'-'Xu-___li__c_u_ 

I    I  itxaju i-i-j'Oj-ii/i  i  <i  i  x^jjctaiaoc  ll   h  J. 

»-  ii  n  <^a-aa.<tJ«i  n  oo  ii  <i  <i —«■»-*  a.  cc  ii  _j_i  n 

3  oooo  i    i    i    I    i   i    i  «.  ii  ii   I  ziw  i    I    i    I   ii  ii  ii  -jujuju. 

Q-  LULU  II    II   II   II   II   II    II    II  LL«C  II  <t  II  Uj  II   II    II   II  r-HCMCOU-UOUJ 

»—  acc_x^LL«LL<x^ooooi_._J'«acu.«a.u.«ioooooooooi— 

~D  a.0,CQCO00(/)CLQ.cDaj  _J  .JO0_LO0_Laji_jO0O0UJUJ_U<_><<O0 

u  j.o.xxxxxxa-aa.a.aaa(£^Koda:^oL<a.aj: 


*♦*  *         *******  *** 


< 

a 

X 

QC 

^v 

a. 

ISJ 

ca 

< 

> 

m 

O 

v» 

co 

1- 

ca 

a. 

a 

* 

1- 

* 

o 

u 

mm 

* 

UJ 

\- 

CM 

a. 

o 

** 

oo 

OO 

X 

\- 

Ul 

H 

co 

oc 

QC 

Z 

Z> 

o 

"^ 

a 

CD 

oo 

s 

t-H 

■^ 

* 

O 

CO 

O 

< 

00 

•-4 

z 

_J 

a 

UJ 

LL 

< 

-L 

o 

u. 

a 

* 

QC 

QC 

LU 

z 

LL 

«^ 

O 

*-»1 

D 

<t 

LL 

■v 

LL 

uc 

O 

►— i 

_l 

<l< 

oc 

CJ 

< 

u> 

oca 

>- 

<[ 

CO 

* 

> 

oo 

»- 

CD 

<coz 

<r 

LL  UJ 

>-i 

oa 

»  2U      CO 

u< 

DO 

> 

a. 

1-X  +  IVJ5T 

X 

** 

•-H 

CQ^ZOCQ 

OO 

1- 

UL 

<TCDU'CDQ- 

LL 

•  • 

►— 1 

o 

VQ.-Q.D 

C 

rlCM 

V) 

***o* 

mm 

Z 

oo 

<<«**cu 

OO 

— *■— » 

Ul 

LU 

ilXUOJU. 

LU 

__Z 

OO 

i> 

Zilaai 

> 

I   ll  i 

I— 1 

>l<ia.cDi 

H-* 

OOl/J 

LU 

1— 

OW  1  «.l_l 

1— 

II  II 

1- 

< 

ii  ii  ii  a  ii 

< 

r-IlM 

_> 

> 

OICOCO  ||  UJ 

^ 

mm 

Q. 

eH 

OZO-rsISC 

•— 1 

-i_C 

£ 

IX 

ajcocnajaj 

OC 

i  H  i 

a 

LU 

a-CLsrxs. 

LU 

OOOO 

o 

O 

OUQOU 

a 

110 


ouoo 

oo 

COCO— < 

CMfO 

ooooa>o,> 

o»o> 

rHw*r-*ft 

.-(— < 

oooo 

oc 

l/)t/)t/)V) 

c/)«/) 

U.tULUMU« 

IWIW 

uiu>trt\/> 

</>VJ 

oo 

oo 

o 

oooo 

i-HCM 

m>o 

oo 

(Mro^irv 

oo 

oo 

o 

■■ll  It  ll'l 

CM  CM 

cm  cm 

CM 

CM  CM  CM  CM 

oo 

oo 

o 

oooo 

WIS) 

WW 

w 

t/7toi//c/> 

UJUU 

IUU1 

u< 

UJ  !__>—.< 

1/1U1 

t/)V) 

«/i 

(SH/)  «/>(/) 

ooo 

o 

o>©-< 

tr* 

•-HCMCM 

C\J 

cm  cm  cm 

CM 

ooo 

o 

www 

V) 

IWUJU.' 

U-< 

to  «/>«/> 

t/J 

CM 


00 


r-l* 

en  • 

2"f\J 


O  < 

CD  K 

a.  ui 

O  X 

— —      M      *  I- 

oo    o    — 

iHCM      CD      <  O 

</>(/)      Q-      I—  I— 
33      O      UJ 

+  +      *      X  I- 


— o 

org 

f«C/5 
t03 
3* 

*xO 

in—) 


-<cm    <r    # 

<_C-      I-      <J- 

oo    iu    -j 

WW        X        < 

Z2     H      + 
«     *      O 

i-»-    m    -J 

"Vs.      _i      < 

-»—  <    «—  w 

«  l     *  to 

i-t-  o    «r  < 

ujui  _ji-:_  z: 

XI  <o<r  < 

Hi-  «--i  l  »s 

*  *      ■*  cu—  — 

m«0      CIO  •     isir>i 
_j_j     _:cocm     mco 

«— •  <>»  +        <_Q_ 

I  +m  I  u-a:     ox 
00<o:u-<nj  +  o 

— I— I  _)<•-<  CC<M  + 
wwQ  CD  03  Q.  X  CQ  w 

— +  a<»«oa.i 
**i>»***  +  ao 

iHCN4u_mct^-r«g#  + 
S_  :__)_:<—:  li-CUrvl 

<<tJ.<tcu<n^<i.a3 

I  luiauu 

II  ii  «-  n  a  ii  a  ii  u 

isIini  ||  ini  ||  rs  ||  rxw 
U.<lNlu.rsl<rvlcO  II 

m  m  co  w  co  w  w  a.  m 

XXXIQ-XIXivl 

UUQUOODQU 


o 

Ul 

a. 
w 

UJ 

C£ 


UJ 

<-> 

cc 

u 

LL 
UJ 

> 
< 

UJ 

X 

u. 
o 

to 

UJ 


< 

h-l 
IX 
UJ 

a 


+  — 

—CM 
i-«C_ 

oca 
a— 

w_: 

<H 


— < 


HUI 
UJX 
XI- 
I-* 

in -J 


CM 

or 
o 


w 

X 

to 

CO 

* 

t— t 

a. 
* 

o 

X 
e_ 

# 

CM 

o< 


CO 
1/5 

< 


I  CJ 


CMUJ 

+  X 

<t  +  xx 


a. 

«_i_Jccaji/)»-i—  x 

_j«co  ens  3 1  x  x  o 
<wcDO-mtoi^i-  \~\-  + 
<-*  co— <oo<«x 

*  >oa#  —*  #  cotoai- 

_j<m_JCM  •  «ooox 

<i-j<*oo  +  +  +  o 
-»*<i**  t  +  iiit 

*  CM*  >T>UUH-I-I-X 

r-CsCn)  —  -*  _J_>LLU.LL.|— 
_<l_.<l  •U.U.CDOC-  CD 
■C  |  <*  I  IM  |  |  J. XXX 
II  II  II  II  |  II  II  <_1  — )_)<_) 
XXXX  ||  XX  II   II    II  — 

_.<_.<»-  U-<|- 1-  I- 1 

co  an  w  w  3  a.  a.  ca  w  a.  |— 

IIIIWXIIIIS 

quqgouqqudu 


CD 


O 
UJ 

a 
to 
u< 

cc 


o 

<x 
o 


> 
< 

UJ 

X 

u. 
c 

UJ 


CD 

co 
a. 
o 

# 

cu 
<I 
I 

II 

CD 


i-H  cQ 

cc  a. 

uj  x 

o  o 


w 
w 

< 


CO 

2. 

< 

to 

X 

o 

COCU  + 

_:_:_i 

CD  CDS 
Q.—.U. 

DOW 

**x 

r-ICMO 

—1—1  + 

<1<I.CQ 

o-iluz; 

W  to  CD 

<.o.u- 

I  I  X 

II  II  u 

coco— 

Z.SL  II 

ll<lO 
ujtoz; 
XX»sl 
OOO 


> 
* 

X 

* 

CO 

m 

<.w 
+  w 

CM< 

ur 

CU4 

*o- 

IMU . 

I    O 

.  II    II 

a.  a. 
c_a 

O-lsJ 

oo 


x 
a. 

o 

cc 

"v. 

to 

UJ 

o 

QC 
U 


> 
< 


u. 
o 

</) 

UJ 


I— I 

cc 

UJ 

o 


< 
I- 

UJ 

X 

I- 
c 


Ul 

a 
c^ 
Ul 

CC 


w 


Ul 


a 

2. 


O 

I- 


LL 
O 

to 


o  < 

II  J* 

X  t- 

a.  cc 

M  UJ 

— i  o 


*^ 
•  _J 

CM< 
w* 
s»CM 

m_i 


<UJ 

*w 

_l  + 

<CM 

ujen 
toz 

I  w 

t- <s 
cnX 
Zh 
t-O 
</!•— • 

COX* 

XI f 

♦  O-J 

lM< 

1-3* 
3*X 

toroto 
O-J -J 

+  <ra  — 

_.*—  X 

XCCU-*  I- 

X     (-<<^«~  < 

0-X<_J_J  •  to 

u-|-C0_._.CM  Q. 

co<a-#w+     q 
Icoo_>*c_     + 

DI  +  -I_<       I 

*uia<_     i- 
ir>*i— -*coco     u_ 

— J  NULL    •auc—wto 
«.— ICUCMO.—    •— - 

I  «!._.   I     I     I   UL > 

II  II  — >  ||    II    II    II  — > 
XX  II  XXXX  II 

»-i-xi-Hr-»-x 

X<lf--IX<Qi|- 

cd  cGma.  ww  a.  w 
a.  a.  a.  a.  a.  o.  a.  a. 

QQUOOUQU 


*** 


##* 


#*■«■        *** 


111 


o 

OCJ 

t_> 

0* 

>3-m 

r- 

cm 

mm 

m 

CM 

cm  cm 

(M 

O 

OO 

O 

00 

COV) 

V) 

UJ 

UJUJ 

u.< 

0O 

oooo 

oo 

oooooo 
CT>o»-*tviro»t 

cm  cm  cm  cm  cm  cm 
oooooo 

000000000000 
UJUJU.IU.UJUJ 
V>0OO'J0'(000O 


o©<j 
o-<cm 

cm  cm  cm 
ooo 
oooooo 

U-'U'UU 

wv>v> 


O 
u 

«l 
I- 
UJ 

X 

>    o 


>s  a. 

—  v. 

X  3 

o:  oo 

Q-  I- 

Q-  Z 

Q  UJ 

+  X 

I  u 

I-  s: 

00 

O-  I 

+  I- 

X  !-• 


a.  u. 

a  o 
o 

+  OO 

X  UJ 

H  > 

co  »-• 

Q.  I- 

tJ  <l 

«-  .* 

II  — 

X  at 

I—  iu 

X  Q 


O 


O 

UJ 

1/5 
UJ 

QC 

"V 
3 

oo 

t- 
z 

UJ 

X 

o 

2. 


u. 
o 


<i       m 

*        z 
a: 

< 

CD 

CD 

a. 
* 
m 

+ 


> 
* 

(VI 

* 

* 
oo> 

*•-• 
t/)<t 

m>» 

(MO- 

^2;  to 

ttO  UJ 

•u.  > 

CMCJ  »-1 

I     II  I- 

III  <. 

U-(-  ^ 

2.0    >-i 
OX      oc 

CLt-       UJ 


o 

X 
co 

3 
CO 


O 
CO 
< 
* 
(X 

< 

CO 
CO 

o. 

I 

M 

a 
co 

rsj     o. 

U-rsjQ 
co<r# 

XCOcO 

ox< 
#CJ~» 
m#* 
—jvoa. 
<».— iu 
l  •w.x. 
ii  ii  ii 

MrslfNI 
U_<tCO 

UJCUQ. 

0.0-0. 
(-jCjO 


-co 

m  • 

ZCM 

*CM 

O-J 

•< 

(MUJ 
>—O0 
>^< 

rH  + 
_JCM 

<«■ 
ujz 

OOt-4 

<O0 

I  "V 

Hi 

ml— 
ZQ 
i  ii  i 

"v#rsia 
kcjco 

l-COOUO- 
O—IU-O 
i— Q-U* 
3~*< 
##U-O0 
U-  —  O0_J 

oo  »_io. 

_J(MQ-# 
Q-  +  #  (M 

wQH-t-J 

♦  <-!<. 

—  CtC0<UJ 

t-<ttlli)UJ 

ii>ajQ.oo<t 
j_aj-»<l  | 
cuu.  I  ii  ii 
3  II  ||  tsirsi 
ajivjrvia-a. 
*  u_<u.< 
•oo  oooooo 

CMO.Q-CLCL 
•QQOQ 


o 


X 


Q. 

> 

> 

00 

> 

a: 

> 

UJ 

•-H 

"V 

•— • 

CL 

< 

3 

< 

^ 

V 

>» 

3 

^» 

00 

*m 

CO 

1- 

«« 

ISl 

1- 

s: 

Z 

• 

U. 

Z 

«i 

UJ 

(n— 

Q. 

UJ 

O0 

X 

+  z 

O 

2. 

a. 

0 

a.^; 

+ 

U 

a 

2. 

o-a 

INI 

2. 

+ 

OOO')      <ico 

CO 

CO 

X 

<tCC       U'_l 

O- 

X 

COCO      2. 

O 

a+     anj 

MM— 

c 

o 

X21C0U. 

H 

—a:    <ct 

O.O.M 

a 

1- 

ecu.  5:00 

•■H 

»^u-.cio:# 

U-«M 

+ 

»-< 

o.oo<a 

a 

• 

<aoi#o 

O0o0O0< 

rvi 

a 

OXOOQ 

z 

I-Ql~CLN< 

a.  a.  a.  co 

CO 

*OI  + 

u. 

> 

ujoo^»**u- 

oooa 

a. 

U- 

o.#  a  co 

a 

0 

±.<icccc#  a 

+  +  +o 

o 

o 

<_>U-#2. 

0 

t—  QCO*)r>i>z 

MMM  + 

+ 

XOO<CO 

00 

a 

-ft-*  a  1  *uj 

U-^U-ivl 

rsl 

UJ 

■«•  -Jl/JU. 

UJ 

acot-iin 

ifMact-4  +  a.Lj* 

v>ioooiL 

OO 

-> 

CUU IU. 

^ 

»-i»o»\ir" 

Hxa  »*X<J 

LLQ.Q.UJI.JQ. 

1—4 

•0.  1  u.o 

1— • 

O 

X    •   "i-ioO    •  1  H-  «-»Qt.cc<_» 

uuaa 

•LJ 

1— 

II  II  II  — 

r- 

• 

O1-*    • 

UlMU    •  11  X** 

II  11  -uo- 

< 

cucocu  11 

< 

O 

oc  11  11  i-«r*  11  —irvio+m  • 

MM  ||    || 

1    II 

-> 

2.2.£UJ 

^ 

II 

UJ'Q.Z'-l 

•Qt<l  II  <I-J    «IM 

a»<MislMM 

I— « 

CQU.-42: 

1— 1 

X 

OUJ<  1 

cm  a.  11  cou.>0  11 

oooococaa-X 

ct 

LLO-JOOX 

a; 

a. 

acta.  11 

11  o)cc_io  11  11  u 

Q.o-a.o.o-1— 

UJ 

ao-o-i— 

UJ 

a. 

3<ooacaKoooziao 

UQQUUQ 

Q 

QQOO 

(2 

a 

0:  at  qcxncc  u  cjcuj^cja 

*** 


*** 


**♦ 


***    #** 


112 


•■HCvl 

occc 

DO 

zz 
<< 


(vjcvi 
** 
** 

-«CM 
-J -I 
CQOQ 
+  + 


r-lCM      t-l 

mm    oc 
o 


CM 

oc 

Q 


O 


o 


D 


o 

»- 

cc 

2 

to 


a 

2. 


c 
o 


coco 


co 


Ul      i-t 


o 

• 

->  o 

a      • 

*  t- 

w  — i     -J    i- 

QO   •    »Z     < 

CC     •_jl-H»H       ^ 

II  l-l>  1 1-    t-t 
a  ii  <-  it  z    oc 

QCZU-XO     LU 

u.a.>-"ao     a 


XI  * 

Q.O.  • 

**  C\J 

(JO  — 

**    ^ 

DO  .-4 

J-X  2 

CX.CC  CD 

**  * 

mm  t-i 

**  —i 

*#  «u 

— —  * 

ooao  • 

CNJCM       CM 

•   •  M«-» 

cmcm<* 

XXOCX 

hhUQ. 
UQt* 


tvj 
CD 
Q. 

X 
Q 
* 


co 

M 

—  X 

OCX 
Qm 
w< 

—      O       ZN 


(M 

CD 
# 


-*2ll^ 

I  I   CC    I 

II  II  u  II 

»4r«a  II  IX 

U-<rviu. 
COCO  CO  CD 

cccccccc 

□QUO 


CM        «-» 

_J    + 

CD       O- 

*  «. 

•  CD 
CM  CD 
wINlO. 
*<* 

xcc< 

Q.U  — 

*  ** 

CM^I-I 

<aia 

I  u~  I 

II  CJ  II 

M  »rM 
<moc 

(DOJQ. 
OCCCaC 


cors! 
I  -J 

2a 

COOC 

*  ♦ 

— rsl 
fMCCCC 
*QCL 
*«~0C 
_JCOO 
LUO  + 
OC<_>M 
>  — CO 

**a. 

iiNO 

L-)l_J  + 
*  P-fO* 
2.U-1" 
U.OCC 

II  ii  a 

rvliNl— 

U_l  II 
>CLIN1 
LL  QCOC 

UQQ 


Ul 

X 


oc 

■v. 
3 

CO 

I- 

iC 

ui 

2. 


C 

a: 

LL 
O 

CO 

> 


^> 

QC 

ui 


—CM 
•"•CC 

oca 
a— 

z< 


Hltl 

mz 

HHCO 

co-^ 

m_j 
_j< 

— * 

ICL 
0-* 
*0 

o* 

*0 
OX 
ia 
oc* 
*m 
m* 
** 
*«~ 
— oox 

COCMI- 
CM    •< 

•  <MCO 
CM^OC 

MD 
Xl-t 
hOI 

CJl-lD- 
H2U- 
2— CO 

«■»  I  oc 

II  II  Ul 
J-X  II 
l-t-X 

cococo 
acococ 

UQU 


_j< 
<* 
*CM 
CM* 
** 
*CM 
i-l_l 
—ICO 
C0  + 
+  — 

—CM 
r-«CC 

oca 

Q  — 

— z 
z< 

I-* 

*    • 

•CM 

CM— 


in. 


<* 

*CM 
H3 
2CD 
CD* 
■R  CM 


X 
X 

H-4 

< 

X 


OJ  Q. 

cu*  oc 
*  »x    o 

•CMH  + 

CM— <  X 

— #  CO  »- 

*>— cc  cc 

HID  OC 

I0_  +  Q 

Q-*i  + 

*<MI-  I 

f-ISCU.  I— 

<i  I  a  'a 
II  t|  ugu 
XI  II  II  — 
l-l-XXII 
LL.<I-I-I 

CDCQOJU-H- 

aCocoCacoC 

QQUQa 


CD 


oc 

2 

CD 

co    s: 

I-  CD 

~<£.  O- 

UJ  Q 

S.  * 

U  h-i 

3L  X 

a. 
-i    * 


x 
o. 


o 
oc 


LL 

o 


CO 


cc 

LU 

a 


a 

_JX 

<x 

l  •-» 
co< 

INI>^ 

— ai 
*». 
aiac 
«*u. 
ii  ac 
ajQ 

2.11 

OCCD 

Q.Z. 

ocac 

DO 


CO 


2: 


c  •  •« 

OC  fMCMt 

LL  XIC 

O  1-1- • 

CO  1-n-iC 

cu  jijex. 


< 

i—i 
oc 

Ul 

u 


rHCM 

22 
CD  CD 


2 

COX 
MX 
*•-< 

ocv. 

o— 
-»I 
za 
i-«_i 

2COQ. 

co  I  oc 

>20 
<-<*co* 
cox>i 

<Q-*0. 

**— oc 

— CMOCLL 
0*OQt 

**i_j*«-»a 

CMfMQ-<_JCO  + 
**<  I  U'OI 
*  *CCCOOCOO- 
rH«MOCM>'-'CD 
_J_IQ<-»*  *  OC 
ffi(U  +  #HIQ 
**IOCQO-  + 

H(\ja<i*Qi 

^JCLLCD2.>U. 

<:<J.ajuJu-LH/) 

I  i  aca.  i  ucc 

II  II  CJ  II    II    II  l_l 

XX  II  xxx«- 

aaiaao.  n 

LL<Q.OCO_JI 

majcDCL^-CLu. 
ococ:q:qcu.qcqc 

UUUUUGU 


*#* 


*** 


*** 


*** 


113 


co 

* 
C/) 


CO 

X 

* 

co 

* 

rn 

M 

m 

o 

* 

o 

J^ 

X 

+ 

a. 

«•» 

2 

>- 

2 

CO 

u. 

to 

IM 

H 

a 

INI 

* 

O 

I 

* 

— 

Q 

^ 

— 

CX 

t-t 

2 

CX 

o 

X 

— 

CO 

o 

■•* 

Q. 

• 

CM 

* 

z 

z 

1— 1 

co 

o 

>v 

^» 

t-t 

<o 

tsl 

H 

a 

CX 
O 

CO 

1 

2 

co 

X 

cx 

< 

** 

2 

CO 

* 

o 

s 

i 

z 

co 

>- 

m 

+ 

2 

CO 

►-» 

> 

* 

ro 

o 

M 

CO 

CM 

* 

M» 

o 

+ 

CX 

CO 

«— 

1 

* 

— 

cocxex 

i 

—to 

o 

CO 

1- 

*■ 

2 

* 

CX 

xox 

—•CO 

exco 

-w 

CO 

Z 

_l 

to 

exx 

o    x 

*-»* 

Ctl/) 

C-I<I 

t- 

zco 

<l 

UJ 

UJ 

>CM 

rvix 

1- 

-J— exx 

-jcocx 

o< 

— z 

o 

•-<< 

z 

z 

a 

*o 

*  -< 

o 

lUflNw 

LUUNX 

H- 

«-z 

z< 

o 

viz 

< 

a 

> 

—  CD 

— < 

a 

CXU*<1 

cxc_>#  x 

C! 

Z<£ 

1—  >» 

(-1 

*< 

»-> 

s 

2. 

* 

CX* 

^■N. 

>- 

>o«~*v 

a. 

>  — H 

O 

►-<>* 

C/}~ 

X 

•"•■v. 

X 

«M» 

• 

QC\i 

■«<  — 

*">— 

*  »><l 

*- 

CO  — 

cc 

*U- 

CX 

V)-» 

CX 

«■» 

-J 

CM 

■—* 

>UQ 

o 

•*  so 

u 

•  — "V.S 

*o 

— >- 

MO 

cx 

_J 

HH* 

lO*«CX 

t~ 

CM«t-0>- 

K 

CMIOO  — 

o 

— >- 

o 

coo 

o 

*Q- 

o 

o 

o 

XO 

C3SU.3 

#r<"i<:o 

♦  x<cx 

K 

cxm2 

H 

x  + 

1- 

cxr02 

H 

■~« 

ex 

Q-al 

OCOCXOCX 

cx 

DOU-D 

CX 

O*  U-O 

ISJOCO 

cX 

*CX 

MQiO 

z 

lfl*QI->C2it 

"V 

CX  |  OCX 

>s 

CXC003 

CX 

s  1 > 

CXM 

ro2 

cx 

s  +  > 

cx 

►-* 

LL 

M2#a.*#  +UIO 

2 

oozo 

2 

OOZCX 

s 

ooo 

s>«» 

OCO 

"N. 

OOO 

s 

co 

o 

*  CO.-H>-U,i-« 

i-«~"t- 

*>UJ  + 

*  +  uicx 

2 

o>>  + 

2  ex 

+  CX 

2 

CXQ.+ 

2 

* 

»Du.oai#o 

to 

HO-O 

CO 

r-ICX— 'O 

sooo 

ococxo 

300 

X 

CO 

*#*G+\aux 

f— 

a>*> 

J— 

OO*  + 

«/> 

•  :d>->- 

co_)coo  + 

CO 

oo-a. 

c/> 

CX 

111 

«/j2.z  +  cm  • 

uua 

z 

*U-CJJfc 

*£. 

*>U*x 

UJ 

Ocxa-u 

UJCXO-CX 

UJ 

CXLLCJ 

Ul 

o 

> 

(\iu.u.OQvMcxa.tj 

uu 

A.LJLJV/) 

LU 

2.U.CJ2 

O 

r-HUC-lCX 

cjixzixo 

CJ 

uqq: 

CJ 

>■ 

e-1 

IM   |     |U|   |i|U   I    CX 

2. 

U.WOLCC. 

2. 

U.UCXCO 

CX 

CU   I   »»,>• 

acc-xtcjcx 

CX 

1  —  *- 

cx 

u. 

f- 

O.C-J   II  l_>  II     1 

II     II   L4 

l_i 

II  II    II  cj 

W 

1  —  II  u- 

u 

Ul  II    HO 

<_l  1  »«*>• 

C_l 

II  II  l_l 

l_J 

o 

< 

ti  -jo>-h  ii  cjq— 

2. 

ooo— 

z. 

II  II  ecu 

a. 

II  CJ  co- 

U- II  >  II  o 

u. 

UD- 

u. 

w 

> 

JL<IQ-U-X->XCL   II 

#•>»>•  II 

exexo— 

co.*-^  II 

CX  1  cx«-» 

U-Q.   II 

II 

M 

ql  ii  oacKaDu 

-1 

Q2QO 

_l 

023  II 

> 

C0020 

>Q  II  2  II 

>■ 

020 

>- 

I 

DC 

>o*-  ii  Qiuixa 

_l 

J-tOZ)*- 

—1 

>«OcXcx. 

<* 

•SCXCO^ 

•«CXULvOCX 

<t 

cxvxx 

< 

CX 

OJ 

a.  ex  llllcXcx  a:  cccc 

O 

u-cxexcx 

o 

U-CXCXCX 

2 

z>>>- 

2»cX> 

2 

>>-> 

2 

> 

a 

uqquduuqu 

cx 

oooo 

CX 

UDQU 

co 

<CJCJQ 

fiQUQQ 

V) 

ooo 

V) 

o 

#* 


#*# 


*#♦ 


*#* 


** 


##* 


**» 


114 


CO 

> 
* 


X 
Q. 
# 

2 

> 

I 

Q 


CM 


Q. 

a 

1 

* 

_j 

2 

2 

< 

</) 

«/) 

1 

> 

> 

«/j 

ft 

0 

w 

H-* 

-j 

*m 

ft 

X 

<CM 

«". 

V) 

—1 

a. 

—  * 

_l 

X 

Ul 

ft 

+  ■* 

111 

* 

cc 

2 

*  <  11 

cc 

_l 

> 

</> 

»—<•—• 

> 

UJ 

ft 

>- 

XX 

* 

a. 

CM 

+ 

a.  a. 

CM 

> 

0 

Q 

** 

0 

* 

ft 

-J 

2? 

ft 

CM 

m 

< 

1/5  CO 

• 

O 

z 

w 

» 

CM 

* 

ft 

w 

+   1 

+ 

• 

1- 

• 

* 

DO 

a: 

CM 

O 

CM 

»» 

-i_J 

t- 

ft 

1 

O 

1 

0 

•«•< 

-J 

0 

1-1 

_J 

H-« 

V) 

1— • 

X 

w»^ 

# 

a 

0 

> 

X 

tvj 

J. 

a. 

—  + 

CM 

J- 

w 

a: 

* 

a. 

* 

a. 

* 

X       *CM 

# 

* 

f-H 

^ 

1/5 

<      2# 

■ft 

0 

ma 

O 

O 

O 

zo 

O 

M 

_l></)# 

> 

1- 

z> 

1- 

<-» 

h- 

*Q- 

H 

1      x 

2"*0~ 

ft 

Di  +  3 

■li- 

QCCM2 

-J2Q- 

^* 

#JUM 

u. 

(X 

XsJ-00 

QL 

ma: 

cc 

X2W 

(X 

>U)X 

n<<ri^fn 

_1*»X    • 

0 

V. 

*zz 

"S. 

azs 

^ 

«  1  z 

^ 

w>Z 

CMM^O-ST 

♦  OOQ-CM 

0 

2 

0*0 

2 

X  +  l/5 

2 

ovia 

2 

#»Q 

XXX^X 

C/)XI«* 

* 

>CM  + 

ftnoz 

Q.ISI  + 

CMX  + 

cm  at  at  2  m  cno 

t/J 

O-it-) 

c/) 

a.xu 

«/> 

0*CJ 

to 

zxx 

^<l«l<i^ 

U\       ITUM*  *l/>*  -R-  X 

•— 

UL  +  > 

r- 

a*  ♦ 

*- 

Cti-iu. 

t- 

**CL 

"luimw     rvi<t  •tr\*-*#cc 

z 

>-i-<LJ 

2 

(XrHtt 

z 

*-Z«-> 

iC 

•    •.* 

— 1 1 — 1 1 1 

'*n-ii-«'T>J-0*U>    •  +  _j_j-k- 

Ul 

LJZCX 

ut 

^■zu 

111 

U  I  cc 

111 

lM>T«/> 

cc.ec<xcc.a. 

WUUI     •    •«J-»*X«JCJ** 

Ul 

i. 

II    II  -£ 

2. 

LJ    |   IX. 

2. 

II  II  *=. 

2. 

II    1  -=. 

uiouuo 

•  •    too    »*»-«tt_l.-ICM   • 

CJ 

UQU 

u 

II    II  ■£ 

U 

QUU 

O 

X  II  O 

h-|— 1    I    |— 

uooii  110  *u.a.<tjcjco 

i. 

*->■  II 

2. 

CLOCLJ 

2. 

U.Q.  II 

2. 

Q-X  II 

^^■^^Z 

II  II    II  £2  II  CM  1  Q_  —  ooi/)  1 

U3Q 

02II 

Q2Q 

2Q-X 

Ul 

i—tt— «»—«•— »t— • 

f-4CM2.<cuu.  11  11  11  11  11  1 

II 

j£ 

«.«/>>- 

2 

(XUJtt- 

2 

OCOOU. 

Jt 

iflxa 

^ 

II   II  II   II   II 

22<_l>Oi-«CMrO>*i-l<MX 

< 

2ZZ 

< 

zzz 

< 

ZZZ 

<t 

zzz 

HH 

i-«cMm^tn 

WWJX^UZ^ZZUUX 

>- 

uaa 

r- 

DQU 

>■ 

QUO 

> 

QUU 

h- 
< 

xx*<;xx 

cc 
ui 

■ftftft                 ftftft                 ft**                 **»                 ftftO** 

115 


UN 

o 


< 

=:        of 

<        x 

■v        * 

— >                *— ■ 

o       o 

IX       Of  Of 

U-     rviu. 

+      *  4- 

-•-■          *•■•..«-» 

o    oo 

Z      OfZ 

<    li-< 

a     +o 

=    -*r> 

o 

of     Oof 

• 

w        Zw 

>o 

-o    z    <z 

X      (-1      Qh 

*        V)       I3</) 

X      -K-       Of  ♦ 

a.    a.    —a. 

CLxO>l"~Z> 

■ 

QXU.XHIL 

+  *#*t/}* 

o 

■4-1  «Q#    • 

• 

XQ.rvJQ.CV  CM 

lf» 

*INlwOf>-*- 

XO  +  CU-  + 

H  +  »0  +  *«0 

Oir\x»o  »x 

XX*X(M# 

rH 

h-*X#wI 

m 

OCQQ.X  i  a. 

o 

+  21>-Q.XZ 

* 

in 

mr«JQ0fXO 

o 

i-t 

xo+om+ 

ii 

• 

♦  +r-  +  <r- 

— 

LO 

o 

CD'TXin*  X 

— o 

-J 

2:x*x— * 

©•H 

UJ 

X*0*t-ICJ 

HX 

t-OI- 

KQ-Q-aif-ta. 

X— » 

O  "O 

-^ 

K- 

uo>-a.xz 

tr\ 

»-«/> 

0»-<0 

X 

Z 

in  +  rvioa*  i_i 

X 

Zl_> 

>o^- 

Q. 

UJ 

xma  +  oof  + 

* 

HU 

•»     •      •» 

•k 

o    2:0 

—                ~*X  +  t-«  +  Cfi-» 

00 

»/■>* 

l-Ot- 

UJ— 

r-4        UJf-H 

_J           «--3cQ#COi-'rOOi-«»-< 

s: 

*CO 

Of  II  Of 

Si> 

»-. 

II   OO   II 

_j       i-i<s:xxxx  +  x^>-< 

co 

O         cox 

3I-3 

l-l  •• 

K- 

a  «<c 

— O— — t-i>CC(~K  «  *  l>#  z>* 

Q.           O 

UiO      X  + 

<_)< 

HX 

O 

0fO_l0f 

r-ofcr>xx:rxxofxxcfXi-< 

a       ix" 

auj     +  — 

>UJ>- 

••    » 

-J 

a. 

XXXX   »  »-C>l-t->|-*Z  +  Qf 

+      •au> 

<ca     — -o 

»o  •• 

«"Hl-« 

Q. 

a    too 

•••■•■   «^l<t  +  C_3MC10fOC3^X 

»-«      <\4<UJ 

Of«t       OH 

>  -> 

•»   •" 

Z 

UIU4     >-,u< 

<^««X— v^H  +  Q  ■*•£_)>  +  •-<  + 

xcoi-tcfo 

*Of        i-4J^ 

►   •   •» 

i-TVI 

w 

*_ HU     UUi 

v*^wJjXHtOl  +  Of  U^^Ct 

*u.*w<» 

_l*        Xw 

xinx 

«— >■» 

-S 

oa.         a. 

— 1—1— l-JU-Ot*  XitXHOXrlM 

rNiUi-imot 

_jsqu,i-»*;<-i 

•-OJ   * 

otu 

Of 

a^.vi    xvj 

OtttOfOf  C>OPM*  X*  ^  T  *  0t"«» 

UtXX*U)UXiWMH 

J-  in 

-S.S 

l_) 

o«* 

u>  <J  <j  u  i— r-  tu  ix  *  u»  *  ex  lj  ■>*.  .£ 

an-t^  ♦  u.-vix  i«i-hjw»^.ui 

O-i-Q-J. 

a.  u. 

if 

viO    »VJ 

1— 1— l-»-^ZZ:irM>-rvlof  >  II  x 

0.[DHHXNXrtl-U#    II  X 

vi-h  >a 

c_>o 

UJ 

Z«J»-    «</5H- 

<£^^*£i-«-i<_ll—  rM>oc.— •*£_!  n 

aa-QX  n  x  ii  x*  *  at-  ii 

mi— tft  ->■ 

•-<— )l_)U^->l_l 

,_,.-<_,_  ||    I)    ||  QQUQ  ||  Q_I5 

II   II  _lf-U  IIQQ->OXofl- 

OZOUUJJ 

-J 

wuz        z 

ii  ii  ii  ii  o— i<u  ii  ii  ii  ii  r-i  M  a< 

COCO  II    ||  DQU  II    II  X  I   -*0 

•»-t  •2.Z_I_J 

_J 

>tXJ£      >^ 

vO  1*-  CO  <J>  f-H-HZ  Q.  M  >  Of  Of  Z  Of  >• 

Q.COQXXQXXO  II    IKQ 

OU- O"-1  □  <  <t  _J  < 

IO 

xx^xxxxxxxoexxxx 

QU-JI-I-  INIQ-O.  XX  >•>>-> 

1- 

uo<o 

•  VJO     too 

ujq:  _j       ^- 

s 

1—  »-ir\i     i-<c\j 

_i«*oci— xu 

*- 

U 

Q.Q.I-ZU.J 

z 

Q — l»-m      t-n> 

ZUJZ»-H<a 

a 

QOQ.O  1  OU  I 

<a:r>ofofcf 

u. 

zi— >*-Jof  «-ia: 

* 

♦  WQ.oa.oa 

1- 

ujco's.a.xoa-x 

116 


LIST  OF  REFERENCES 


1.  Naval  Postgraduate  School,  Progress  Report  to  SESPO, 
Pressure  Ratio  Effect  on  the  Heave  Motion  Characteristics 
and  Pressure  Dynamics  of  XR-3  for  Step  Weight  Transients, 
Alex  Gerba,  Jr.,  and  G.J.  Thaler,  January  1977. 

2.  Barnes,  L.W.,  Development  and  Time  Domain  Validation  of 
a  Low-Order,  Constant  Speed,  Pitch-Heave  Model  for  the 
XR-3  SES,  M.S.  Thesis,  Naval  Postgraduate  School,  March 
1977. 

3.  Riedel,  R.G.,  Sensitivity  Study  of  the  XR-3  Loads  and 
Motions  Computer  Porgram  Sidewall  Parameters  and  Forces 
on  Roll  Behavior  in  Calm  Sea  and  a  Comparison  to  Test- 
Craft  Turn  Maneuver  Behavior,  M.S.  Thesis,  Naval 
Postgraduate  School,  1977. 

4.  Department  of  the  Navy,  Surface  Effect  Ship  Aero/ 
Hydrodynamics  Technology  Design  Manual,  SES  Project 
PMS-304,  Vol.  II,  February  1976. 

5.  Department  of  the  Navy,  Surface  Effect  Ship  Aero/ 
Hydrodynamics  Technology  Design  Manual,  SES  Project 
PMS-304,  Vol.  Ill,  February  1976. 

6.  Menzel,  R.F.,  Study  of  the  Roll  and  Pitch  Transients 
in  Calm  Water  Using  the  Simulated  Performance  of  the 
XR-3  Loads  and  Motions  Program,  M.S.  Thesis,  Naval 
Postgraduate  School,  1975. 

7.  Leo,  D.G.,  and  Boncal,  R. ,  XR-3  Surface  Effect  Ship 
Test  Craft;   A  Mathematical  Model  and  Simulation 
Program  with  Verification,  M.S.  Thesis,  Naval  Post- 
graduate School,  December  19  73. 


117 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  014  2 
Naval  Postgraduate  School 
Monterey,.  California   93940 

3.  Department  Chairman,  Code  62 
Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  California   93940 

4.  Professor  George  J.  Thaler,  Code  62  Tr 
Department  of  Electrical  Engineering 
Naval  Postgraduate  School 

Monterey,  California   9  3940 

5.  Professor  Alex  Gerba,  Jr.,  Code  62Gz 
Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  California   93940 

6.  Professor  D.M.  Layton,  Code  67Ln 
Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California   93940 

7.  Mr.  H.  Robertson 
PMS  304,  Code  141 

Surface  Effect  Ship  Project  Office 

P.O.  Box  34401 

Bethesda,  Maryland   20034 

8.  Deniz  Kuwetleri  Komutanligi 
Bakanliklar- ANKARA  /  TURKEY 

9.  Deniz  Harp  Okulu  Komutanligi 
Heybeliada-ISTANBUL  /  TURKEY 

10.  Istanbul  Teknik  Universitesi 
Gumussuyu-ISTANBUL  /  TURKEY 

11.  Fuat  Ozanturk 

Esentepe  Emekli  Sb.  Sitesi 

18.Blok  D.5  Mecidiyekoy- ISTANBUL  /  TURKEY 


118 


cur"  maneuvers. 


Thesis  186344 

0986   Ozanturk 

c.l        Development  of  the 

five  degree-of-freedom 
linear  model  for  the  XR- 
3  surface  effect  ship 
and  investigation  of  the 
roll  behavior  of  the 

craft  in  turn  maneuvers. 


Ihes0986 

Development  of  the  five  degree-of-freedo 


3  2768  001  00091  2 

DUDLEY  KNOX  LIBRARY 


